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CONTRIBUTION TO THE LINKAGE THEORY OF 
AUTOPOLYPLOIDS: I 
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_ A theory of linkage of autopolyploids is developed under considera- 
tion of m loci and * alleles. The simplifying assumption of chromosome 
segregation, which may be considered as an approximation to the more 
adequate theory of chromatid segregation, is made throughout. Random 
mating and distinct, non-overlapping generations are assumed. Under 
these assumptions the problem is determined by three basic probability 
distributions—the distributions of genotypes and of gametes, and the 
segregation distribution. The segregation distribution is assumed to be 
the same for males and for females. The aim of the paper is to estab- 
lish recurrence formulas (which allow to find the distributions of 
gametes and of genotypes from generation to generation, if the dis- 
tribution of genotypes for an initial generation is known) and to inves- 
tigate the limit behavior of these distributions as the number of genera- 
tions increases indefinitely, 

In the present paper (hereafter referred to as I) the problem is ex- 
plained, and the three characteristic distributions are introduced for the 
general case of a 2s—ploid, m loci, and 7 alleles. Recurrence relations 
are established for tetraploids, s—2 and m=—2 loci, while the recurrence 
relations for the general case as well as the limit theorems will be given 
in the second part of this paper (hereafter referred to as II). 


Introduction. 

A linkage theory for tetraploids, 2s—4, based on “chromosome 
segregation,” has been presented by D. De Winton and J. B. 8S. Hal- 
dane (1931) for two linked factors, m=2, and r alleles. The case, 
s=2 and m=3, is also briefly considered. The same problem has been 
studied by F. W. Sansome (1933) under the assumption of “chroma- 
tid segregation.” In a more recent paper R. A. Fisher (1947) pre- 
sents a linkage theory for tetraploids and hexaploids also under the 
’ assumption of chromatid segregation. The approach is similar to 
that in K. Mather’s (1936) study on segregation and linkage in 
autotetraploids. hs 

This author has previously considered the linkage of normal, 
that is, diploid organisms, s=1, with m loci and r alleles (Geiringer, 
1944; 1945). In these papers, as well as in the present paper and in 
the papers just quoted; random breeding and distinct, non-over- 
lapping: generations. are considered. By means of an appropriate 
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“linkage distribution” general results have been obtained. Recwr- 
rence formulas have been derived which enable us to compute suc- 
cessively the “distribution of genotypes’ as well as the “distribution 
of gametes” from generation to generation. The main result consists 
of a general limit theorem as n, the number of generations, tends 
toward infinity. 

On the other: hand, following up an important paper by J. B. 8. 
Haldane (1930) the author presented (Geiringer, 1948a) a mathe- 
matical theory of autopolyploids in the case of 2s chromosomes, one 
locus, and 7 alleles. A “segregation distribution” is defined which 
contains random chromosome segregation as considered by J. B. S. 
Haldane as a particular case. Again recurrence formulas are estab- 
lished and integrated, and the equilibrium status which presents it- 
self, in general, as 2 > © only, is determined. A further paper which 
will be published in the November issue of Genetics deals with chro- 
matid segregation of tetraploids, hexaploids and octoploids form=1. 
The segregation distribution which will be introduced includes, as 
particular cases, the segregation distribution, (s.d.), corresponding 
to chromosome segregation (Geiringer, 1948a) and random chroma- 
tid segregation as assumed by J. B. S. Haldane (1930), and R. A. 
Fisher’s and K. Mather’s (1943) segregation theory. Consequently, 
the results obtained in this paper—recurrence formulas, integration, 
and limit theorems—apply to all these situations. 

In the present paper, returning to the assumption of chromo- 
some segregation, we consider the general problem of m loci for a 
2s-ploid with r alleles. An adequate setup must include as particu- 
lar cases the problems of diploids of m loci (Geiringer, 1944), s=1, 
general m , and of a 2s-ploid and one locus (Geiringer, 1948a), m=1, 
general s. Moreover the case of a 2s,-ploid with m, loci, (s, = s, 
m, =m), must appear as a “marginal” case of the more general case. 
More explicitly, if by an appropriate marginal summation, applied to 
any of our distributions, for example, to the distribution of gametes 
or to the s.d., we reduce the m to m, and the s to s,, the marginal 
distribution so obtained must be a distribution of gametes or a s.d. 
for a 2s,-ploid with m, loci. 

The theory is based on the definition of a general s.d. which is - 
given in section 2 for m=s—2, and in section 4 for general m and s. 
This s.d., which, of course fulfills the above requirements, contains as 
particular cases random chromosome segregation, as considered brief- 
ly by J. B. S. Haldane (1930) and De Winton and Haldane’s segre- 
gation theory( see section 3 of the present paper). The main result 
is a surprisingly simple and general limit theorem which also com- 
pletely settles the question of the existence of an equilibrium status 
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(sections 3 and 5 of II). It is based on recurrence relations (section 
1 of II) between the distributions of gametes in two successive gen- 
erations. (From the distribution of gametes we easily deduce the dis- 
tribution of genotypes, if random mating is assumed.) The recur- 
rence formulas form a system of quadratic difference equations whose 
constant coefficients are the values of our s.d. 

As has been said, the theory applies to chromosome segregation 
only. It seemed necessary to present the theory under this simplify- 
ing assumption before dealing with the very considerable additional 
complications of chromatid segregation, considered, so far, by the au- 
thor for m=1 only (Geiringer, 1949). Chromosome segregation forms 
an approximation to the more adequate assumption of chromatid seg- 
regation. It is also assumed that the same s.d. applies to males and to 
females. The desirable generalization to different s.d.’s for the two 
sexes may not be too difficult (see Geiringer, 1948b where this is 
carried out for s=1 and general m). It is however not attempted in 
the present paper. The ideas developed in this paper may seem too 
abstract and formal to the majority of geneticists. The difficulties of 
some problems of modern heredity are quite considerable and seem to 
call for a mathematical approach. Since in science there is a mutual 
interaction between observation and theory, it is hoped that such an 
approach may, sooner or later, prove useful to theoretical genetics, 
although it emphasizes only one side of a complex problem. 

Let us finally point out that the linkage theory presented here 
forms the framework of. a theory rather than a specific theory, inas- 
much as no particular assumptions are advanced regarding the para- 
meters which enter into the definition of the s.d. The theory is to be 
completed by biologically adequate “models” of a s.d. such as in R. 
A. Fisher and K. Mather (1943), or in H. Geiringer (1944). In this 
sense Haldane’s linkage theory constitutes a particularly simple 
model. The results of the present paper apply to a great variety of 
models. Thus they describe, e.g., the limit behavior of the basic gene- 
tic distributions for De Winton and Haldane’s linkage theory. 


1. Notations and Definitions. 


As in the introduction, let 2s denote the number of chromosomes 
in the zygote; m, the number of loci; 7, that of alleles. We assume 
that the concepts for m=1 are more or less known (see Haldane, 
1930; Fisher and Mather, 1943; and Geiringer, 1948a). We start 
with the consideration of a tetraploid with two loci and two alleles, 
m=s=r=2. Our first task is to enumerate the possible gametic and 
zygotic types and to introduce a suitable notation. It we attempt this 
we see immediately that there are various possibilities of differentia- 
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tion. We shall mention four of them, beginning with the considera- 
tion of zygotes. 


1) One may assume that a zygotic type is satisfactorily char- 
acterized if the four values are stated which correspond to each locus, 
such as (AAAaBBBB) or (A‘aB*). This is, in fact, sufficient if, as in 
J. B. S. Haldane (1930), no linkage is considered. 


2) De Winton and Haldane, in their paper on linkage, show in 
their notation the four chromosomes which constitute a zygotic type, 
denoting such a type, for example, by (A4B,AB,Ab,ab) where the 
chromosomes are separated by commas, and where the order of the 
chromosomes plays no role. In this conception (AB,Ab,Ab,aB) and 
(AB,AB,Ab,ab) represent different types, while in notation 1) they 
both appear as (A*aB?b?). 


3) We may choose to show the maternal and paternal heritage 
of the zygote in our notation. The reason for this will be discussed 
at the end of this section. This may be done without indicating the 
four chromosomes. We may simply write the “maternal values” 
before a semicolon and the paternal ones after, for example, 
(A2b?;AaB?). If we denote the whole maternal heritage of an or- 
ganism by x, the whole paternal one by y, we shall assume that 
(x5y) = (y3x). 


4) Finally, parental heritage as well as chromosomes may be 
shown. A type may then be denoted by (AB,Ab;AB,aB) and we 
assume, as before, that («;y) = (y;x) and that the order of the two 
chromosomes before or after the semicolon plays no role. 


Let us count the number of different types under these various 

assumptions. Without any “symmetry relations” there would be 

*&—=256 different types. According to the four assumptions this num- 
ber reduces considerably. There are under the 


first assumption n,—25 different types; 
second ” Ne—o0 é SERV 
third sit N3—45 a ies 
fourth \” NAH=0 fe x 


These numbers are easily understood. We find all possible types in 
the first case by symbolical multiplication of (At + 443a + 6A2a2 + 
4Aa’ + at) by (Bt + ---- + b*). Hence there are 25 different types. 
The coefficients, for example, 24(A°aB?b?), indicate how many of the 
original 256 types are lumped together. In the second case we shall 
denote, for the moment, the four chromosomes by u, v, w, and s. 
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The possible types correspond to the terms appearing in the ex- 
pansion of (ututwts)*, hence comprise 35 types, and the coeffici- 
ent of each type appears also in this polynomial formula. In a similar 
Way we may establish the number of types for the third and fourth 
cases. 

The enumeration -of the gametic¢ types is simpler.. Here no ana- 
logue to notations 3) or 4) exists. We merely may or may not show 
the two chromosomes in a gamete. 


1) If they are not shown, we assume that the gamete is 
completely specified by stating the four values, such as (AABB) = 
(A?B?). There are then obviously nine different types. 


2) If chromosomes are shown then the types (AB,ab) and 
(Ab,aB) are distinct and there are ten types. 


Now consider a 2s-ploid with m loci and r alleles. A zygotic 

type is then denoted, e.g., under the fourth assumption: 

(7 oes Ly) A » aH 98) ee 6 Py) : 

yy fetate Ym™ yee Yy °°) Bis Ym"). 
Here x; is the maternal value corresponding to the 7th locus in the 
jth chromosome, and y;%) means the same for paternal values, 
((=1,----m,j=1,--:: 8). There are 2sm places, each to be filled 
with one of 7 numbers. Consequently there would be N’=r?*" types 
if there were no symmetry relations at all. Corresponding to the four 
ways of differentiation we get the following numbers without diffi- 
culty: 


acess il (ns seaitaaa) 
ta pita ( (eat at 28 : 
ot Ce an 
This reduces for s=1 and any m and r to eek which was de- 
noted by Nz in H. Geiringer (1944) ; while for m=1, and any s and 


(Mig APSE | ‘ ae 
r, we find the number N,’ =| ae ), as in H. Geiringer 
(1948a). | 


2) This leads to 


“ee (Bain sien) 
ie 28 f 


which for s=1 gives N, = 4 ym (7 + 1) as in H. Geiringer oa : 
while for m=1 we again find N,’. 
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aC ST ad 


s-r—1 
For s = 1 we again have N, and for m=1, with R= ( : ) , we 


have N,' = 4 R(R+1) as in H. Geiringer (1948a). 
4) Finally 


ne = pL oe 
eae" +8 *) {C is Ei 
Ss s 


which for m = 1, any s and 7, reduces to N,’ and for s = 1, any m 
and 7, to N;. 


If 7; denotes the number of possible alleles at the ith locus, 


(i=1,---- m), the corresponding expressions are: 
per) v; oe PASS —_ 1 
ee: ( ) 
4 —— 1 
1 {gy 0 ODA En ar 2s—1 
2) ( ; 
2s 
1.... s ao Ti eae 1 Ss + i aime 1 ye 
3 6OC«@it ( ) takes the place of | ; 
a Hie ho il r—l1 


A) 1, 1, -+++ %» takes the place of 7”. 


We may however assume without loss of generality that there are r 
possible alleles at each locus. In fact, we shall denote by 7 the great- 
est of the m numbers, 7;. To assume this same r at each locus 
amounts merely to the assumption that the genetic probabilities of 
certain types within the distribution of genotypes are zero. It is an 
unnecessary complication to consider, as some authors do, all pos- 
sible combinations connected with different 7; at the ith locus. 

The differentiation between the various chromosomes as in 2) 
and 4) exists neither for s = 1, general m, nor for m = 1, general 
s. In these situations only assumptions 1) and 3) come under con- 
sideration, and assumption 3) was considered by H. Geiringer (1944 
and 1948a) ; (see end of this section). 

Next let us introduce the probability distribution of gametes 
(d.ga) and the probability distribution of zygotes (d.z), or of geno- 
types. We consider again the case m = s = r = 2, the general dis- 
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tributions (which are completely analogous to those form=s=r=2) 
being needed in section 1 of II only. We shall denote by 


2 2 ot 
#4 a Ze 0 al Ohya 


the probability of a gamete of type (A*a?-"B%h2-) in the nth gen- 
eration, where 


So ( i ) ( : Jom (Ang? Babes) —1, 
(ier ae ie (1) 
or, more explicitly, 
p™ (A2B?) + 2p” (A?Bb) + p™ (A2b?) 
+ 2p™ (AaB?) + 4p. (AaBb) +----+ p™ (a2b?) = 1. 


We shall also denote, according to assumption 3), by 


2 2 D) 2 
( ( ) ( ) ) wm ( A%Q?-*1 Bah 2-22 ; Ang2-1 By b2-%2 ) 
U3 Hh 2 Yn Yo 


the probability in the nth generation of a zygote of the type indicated 
within the parenthesis. Here 


Oieeet Osaece Os cce S( =) ( 2 ) ( "| ( i) 
= “ 2 2 X ; Xo Y1 Yeo 
-W™ (A%02-1B%h2-2, Ang?24Bwrh?) =1, (n=0,1,--:-). (2) 
As before, we shall sometimes denote the whole expression 
A” q?-"1 B* b2-*2 

which in w stands before the semicolon, by x, and similarly we shall 
denote by y the analogous expression after the semicolon. The fact 
that we consider (x;y) and (y;x) as identical types corresponds to 


assuming 
wm (x;y) =w™ (y;x), (n=0,1,----). (3) 


Under random breeding a new individual of the (n+1)st gen- 
eration is formed by the random fusion of two gametes of the nth 
generation. This is expressed in the basic formula, 

wn) (x 3y) =p™ (x) p™ (y) 5 (n=0,1,-:--). (4) 


Let us conclude this section by commenting upon notation 3) 
adopted by us, which amounts to introducing a differentiation be- 
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tween the paternal and maternal heritage of an organism. This dif- 
ferentiation is assumed merely with respect to the formation of 
gametes. In the simplest case, s=1 , m=2, we know that there is a 
difference between an organism (AB;ab) and (Ab;aB) with respect 
to the formation of gametes, i.e., with respect to linkage. The former 
will transmit (Ab) or (@B) with probability ‘ce’ (recombination 
value) ; while the latter will transmit them with probability (1—c). 
Thus if linkage is considered there are not nine but ten types for 
s=1, m=r=2. For larger “m’” this differentiation seems unavoidable 
in connection with linkage (see Geiringer, 1944). It seems to the 
author that, if s > 1, the theoretical consideration of linkage calls 
even more for this differentiation.} 

This differentiation is indispensable in order to find and to un- 
derstand certain regularities whose proof forms the subject of the 
author’s genetical investigations. It is, in fact, possible to show that 
in rather general situations certain biological laws hold which are of 
a clear, simple and uniform structure. They are concerned with: 1) 
recurrence formulas which show how the basic distributions of 
gametes and of zygotes in the (n+1)st generation are determined by 
the distributions in the mth generation; 2) the characterization of a 
possible equilibrium status, a problem that includes the study of the 
limit behavior of the distributions as n tends toward infinity. For 
these purposes, which are of a biological nature, the more general 
notation 3)—far from introducing an unnecessary complication—is 
of decisive importance. As a very simple illustration of this consider 
equation (4). If we do not show paternal and maternal heritage we 
would, for example, denote by 


(4): (S) anew 


the probability of a zygote of type (A4*aB2b?) and we would have 


240) (A2aB?b?) = 2[2w™ (A2B?; Aab?) 
+ 8w™ (A?Bb;AaBb) + 2w™ (A2b?;Aab?) ], 


and therefore, corresponding to (4), 


Bui" (A5aB?b?) = p™ (A2B?) p™ (AaB?) 
+ 4p™ (A?Bb) p™ (AaBb) + p™ (A*b?) p™ (AaB?), 


instead of the clear formula (4), thus 


_ _ tIn Haldane’s linkage theory (section 8 of the present paper lyploid 
linkage is reduced to the consideration of diploid none aE ce ce 
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aw (mt) (A2B2 ;Aab?) — p™ (A2B?) p™ (Aab?) ¥ 
w) (A?Bb ;AaBb) =p (A2Bb) p™ (AaBb), ete. 


For greater values of m and s these complications would become still 
more marked. 

Another remark regarding the differentiation contained in nota- 
tion 3) is the following: Our final recurrence formula (section 1 of II) 
is for gametes, which means we compute p'"*) from p™, hence no- 
tation 3) in no way enters into this final result. A similar remark 
applies to our limit theorem. On the other hand, if we want a limit 
result or a recurrent computation for the d.z. this follows immedi- 
ately from the respective result for the d.ga. by means of (4). Final- 
ly, it is extremely easy to eliminate this differentiation 3) after it 
has served its purpose if we are concerned with a problem in which 
it seems irrelevant. We have merely to omit the semicolon and lump 
together the values separated by it. 


2. The Segregation Distribution for m=s—2. 


Let us again start with the consideration of a tetraploid with 
two factors, the case mainly considered in the literature of our prob- 
lem. The fundamental definition is that of the s.d. Our definition 
will again be based on assumption 3) where we show the parental 
heritage but not necessarily the various chromosomes. The number 
of alleles may be arbitrary. Consider a zygote of type (@,a.bib2; 
a,0,b:b,) where a; denotes an allele with respect to the first factor, 
b; one with respect to the second factor. We introduce 


( aivell 2 yy A(&1 €2) } (c,=03;1,2;e=0,1,2), 
&1 Eo 


the probability that «, maternal a-values and e, maternal b-values be 
transmitted, (and, consequently, 2—e«, paternal a-values and 2—é., 
paternal b-values). Thus 4(00) is the probability that no maternal 
value will be transmitted, and, consequently, the four paternal val- 
ues, d;, ds, b;, bs, will be transmitted. In the same way 4)(10) is the 
probability that either (@,a3;b3b.), (d:@sbDsbs), (Q2@sbsb4), or (a20,bsb4) 
will be transmitted and we assume that these four probabilities are 
equal. Moreover, as in previous papers, we have to assume “sym- 
metry” with respect to paternal and maternal heritage. This amounts 
to puttin 

; : A (er 22) =A(2 — t1 5282). (5) 


Hence there are five parameters with sum one. In fact, because of 
(5), the equation 
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A(00) + 42(10) + 44(01) + 164(11) + 2(02) 


+ 4(20) + 44(12) + 44(21) +2(22) =1 a 
reduces to 
21(00) + 84(10) + 84(01) + 164(11) + 22(02) =1. (7) 
It will be convenient to introduce also 
ueae)=(— \ * )ateed, (8) 
By & 


the probability that a fixed set of e, maternal a-values and a fixed 
set of e«, maternal b-values will be transmitted. If, for example, the 
fixed set consists of a, and b,, 


1(11) = Prob [(a:b,a3b3) + (@,b,a3b.) + (@,b,a,b3) 
+ (a,bia.b,4) ] = 41 (11), 


where ‘Prob [----]” means “probability to transmit the expression 
in brackets.” 


We then have 


FFCV (eo 882 


=1. (9) 


Let us now introduce certain “marginal distributions” of our 
s.d. With the above notation for “Prob [.---]” where we agree that 


the order of the four values within the parentheses is irrelevant we 
consider 


Prob [ (a142b,b2) + (a,@2b,b3) + (a1a2b,b,) + (@,a2b2b3) 
I- (a,02b5b,) mi (,Q.b304 )] 5 


This is obviously the probability that, no matter what happens to the 
b-factors, the maternal a-factors, a, and a,, are transmitted. Hence 
it is a marginal probability and equal to 


A(22) + 44(21) + 4(20) =1(22) + 21(11) + 1(20) =4,(2). 


Zz, 
In general ( * Jace) is the probability that, as far as the ith fac- 


tor is ei kited é maternal values are transmitted. Thus for 
e=0,1, 2 we have 
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L(e) =l(e, 0) tealle il yet £ (ss 2) 5 (10) 
and 


iL (e) =1(0, e) + 21(1, «) ani2n e). (11) 
with the relations 
(HO) St(2)5 61;(0)-F 21,(1).+-1,(2) =1; G=1,2). (12) 
We note that each of these two marginal distributions has the same 
properties as the s.d. for m=1 in H. Geiringer (1948a). 
On the other hand consider the following marginal prob- 


ability for the event that a, and b, are transmitted no matter which 
are the other transmitted values: 


Prob [(a1b:d2b2) + (@,b1@2b3) + (a1b,42b,) + (a,b,a3b2) 
+ (a,b,a3,b3) + (a,b,a3b,) a (a,b,4,b) 
We denote this probability by 1 (11). This notation indicates, 
with respect to the superscripts, that the “s” of each of the two loci 
has been reduced by ove. If there is no possibility of a misunderstand- 
ing we may also denote this distribution by lw.(e, 2), (e, = 0, 1), 
since it is a linkage distribution with all the properties of the linkage 
distribution for diploids as in H. Geiringer (1944). We then have 
be ed (22) 2A C21) i 2A (12) 4. (11) 
=ADRPu abet 1(12).i4- (22); 
BLO) 14106 20) 4 (11) + C1): (13) 
bo (01) = 1(O1)_+ 111) +1002) +7612); 
L.(00) =1(00) i+ 7(10) +2(01) + 2(11) ; 
with the symmetry relations following from (5) 
Le(00) =h2(11), b2(01) =he2(10). (14) 
From (9) we conclude 
Lo (00) + L2(01) + he(10) + L211) 
= 21,.(00) + 21.(01) =1. 


From this linkage distribution we derive, in particular, the recom- 
bination probability (r.p.) 
We have thus introduced a s.d. for m=s—=2, which contains the 


five parameters 1(00), 1(10), (01), 1(20), and J(11) which, because 
of (9), reduce to four. The following relations hold for these para- 


meters 


(15) 
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1(00) + 21(01) + 1(02) =1,(0) ; (16) 
1(00) + 21(10) + 1(20) =1,(0) ; (17) 
1(00) '+ 2(10) + 2(01) + 1(11) = he (00) ; (18) 


where 1;(0) is the parameter of the s.d. for s=2 with respect to the 
ith character and ¢,. = 1 — 2l2(00) is the r.p. between the two loci 
as s=1. If we know 1,(0), £,(0), and ¢. from observations on mono- 
factorial tetraploids or bifactorial diploids of the same plant, there 
are four relations between the five parameters. 

Let us finally characterize the important particular case of “in- 
dependence” for our s.d. We call a s.d., [(é: &9), “independent” if 


I (€ &) = L,(e1) L.(e2), (19) 


where, as we easily see, 1;(e:) are the marginal distributions intro- 
duced in (10) and (11). In this case the matrix 


{I (1 &2)}, (ec: =0,1,2;2—0,1,2), (20) 


is of rank one, and vice versa, as can be seen without difficulty by 
means of (5). It also follows that in this case 


le (00) = ke(10) = he(01) =42(11) =F, G2= 9; (21) 


“a 


hence there is no linkage, in the sense of diploid linkage. A particu- 
lar case of independence presents itself if all 


1 
A ey ea)? gad US yy se é—=— 07172 
(22) 


1 
L, =1, ee 1 = as 
(0) (0) 7 21,(1) = 21,(1) ; 


Hence, if the s.d. of the tetraploid is independent, and, in particular, 
af it is uniform, 7.¢€., if each A(e, &2) = 1/86, there is no diploid link- 
age. If there is no diploid linkage the s.d. is independent. 

The definition of our s.d., as given in this section, is again based 
on the differentiation 3). We shall comment upon this after the role 
of the s.d. in our recurrence formulas (section 1 of II) has been es- 
tablished. Let us finally point out that under the assumption of chro- 
mosome segregation a s.d. for tetraploids with two loci cannot ac- 


4 


4 
count for more than ( ) : ( 
Z 2 


) = 36 gametic outputs. If chromatid 
8 


> ie 28? possible gametic 


8 
segregation is considered , there are( 5 ) : ( 
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outputs.; Our aim in this paper is to carry out the consequences of 
a general chromosome segregation theory. 


3. De Winton and Haldane’s Linkage Theory for Tetraploids. 


The linkage theory proposed by De Winton and Haldane with 
respect to the tetraploid primula, m=s=2, is as follows. A zygote 
consists of four chromosomes which are denoted by A,B, , A>B., A;B; : 
A,B, or by K,, K., K;, K,. Starting from observations the authors 
arrive at a linkage theory by the following considerations. In one 
third of all cases K, will pair with K, and, on the other hand, K, with 
K,. In one third of the cases K, pairs with K; and K, with K,, and 
in one third of the cases K, pairs with K, and K, with K;. Next it 
is assumed that for each such pairing ordinary diploid linkage pre- 
vails. That amounts to the following: 


As K, pairs with kK. , 


4 , ; A q mace, 1 — Cx 
A,B, will be transmitted with probability = = —_ = 3 ; 
+? 39 a9 p — Cr2 e 
ASB es 5 ae ; 
3? 3) a9 ob) +”) p ‘a 
A,B, 2 ’ 
33 3) +” BB) 3) q . 
A.B, 2 ’ 


and as K; pairs with K,, , 
A;B; will be transmitted with probability q/2; 


A,B; Joey ” a” oe) /2 : 
ALB hd te ” ” ” p/2 : 
Ali: 29799 ”? 2” ” q/2 : 


Next, assuming independent combination of probabilities; it follows 
that: 


2 


eae @ 
A,B,AsB; will be transmitted with probability Pine 


: t be ob- 
it 28 . 28 gametic outputs there are many which canno : 
egies picasa ea eeceston There are 100 different capa Sei 
Besides the 36 which appear under chromosome segregation, t mee ie ee 
type A,2B;2 (G=1,---' 457 =1,---: 4), and 48 of types A,°B; p OF DA ALE 
2B; ; 


where j <k. 
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4 5 , 1 pa. 
ASE ALB fae ; mr 

, ” ’” 1 Pq . 
A,B,A.B; OP LB , Bape 
Dp Mk So Ne ie eere PR LS, Se 4 ) ii 
A.B,A.PB, ay 3) x) 33 3 ee es A 


Thus, in an obvious notation we see that if 
1) (A,B, ) A.B) x (A,B; , A,B,), 


the four combinations 
A,B,A;B3 
A 1B,A.B, 
A,B2A;B; 
ALB ALB, 


9 
2 


are transmitted each with probability se while 


A,B,A;B, 
A,B,A,B; 
A.B,A;B, 
A,B,A,B; 


2 
are transmitted each with probability ae and the eight combinations 


A,B,A;B, A,B.A;B; 
A,B,A.B; A,B,A.B, 
A,B.A;:B, A.B,A3B; 
A.B.A He A,B,A sD, 
are each transmitted with probability pq/12. Next if 
2) (A,B,,A;B;) X (A,B, AsB,), 
the four combinations 
A,B A.B. 
A,B,A.B, 
A. B2ABs 
A;B;3A,B, 


are transmitted each with probability a while 
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A,B;A.B, 
A,B;A,B, 
A;B,A.B, 
A;B,A,B, 


are transmitted each with probability o and 


A,B,A,B,, A,B,A,B,, A;B;A,B,, A;B;A.B, 
A,B;A2B, ’ A,B;A,B, , A;B,A;B, ’ A;B,A.B, 


all with probability pqg/12. The results of the pairing 
3) (A,B,, A,B.) X (ABs, AsBs) 


can be worked out in a similar way. 


In order to summarize these 48 results we briefly introduce, just 
for this purpose, the following probability distribution. Denote by 
lo. each of the six probabilities that the four transmitted values come 
from only two chromosomes. Thus, e.g., 2 is the probability for 
transmitting A,B,A.B.. In the same way 6l,::, denotes the prob- 
ability that the four transmitted values are from four different chro- 
mosomes, e.g., A,B,A;B, and 24l.,, is the probability that three dif- 
ferent chromosomes take part, e.g., A41B,A.B;. We then have 


Ge2:c-Onm 1 24, — 1. (23) 
De Winton and Haldane’s theory then amounts to assuming that 


loo = ge) bin =—, bua =—, (24) 


in accordance with (23). 

Let us indicate a few further consequences of this theory. Our 
postulate regarding marginal summations (see the introduction) is 
satisfied. We get, in fact, 

1 4 
1,(0) =2,(0) et 21,(1) = 21, (1) ae 
(25) 
{ p a q 
Ls (00) = bye (11) meets lL» (10) = lb» (01) Sob 


Next, it is easy to compute our values / (00), 1(10), ---- 2(02) - (sec- 
tion 2) for Haldane’s theory. We find: 1(00) = L(22y = 2/6, 
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2 
21(01) = 21(10) =21(21) = 21(12) = pq/3,41(11) ae (pr Po as 
1(20) = 1(02) = p?/6. 


The authors point out that in this theory it is possible to dis- 
tinguish “between 50% crossing over and the absence of linkage.” 
That means, in our notation, that for p=q=—1/2 we do not obtain 36 
equal segregation probabilities as we would for random chromosome 
segregation. In fact, if p=—q=—1/2 we obtain 

1 1 
bre hn og? aa Ae 
This is a consequence of the assumption, suggested by the material 
observed, that pairing occurs in the above indicated manner. More- 
over, a uniform s.d., or random chromosome segregation, can never 
be obtained in this theory no matter how we choose the parameters. 
In fact, putting 


(26) 


1 
bos = bu = han = 36 
leads to contradiction with (24). (If q?/6 = »?/6 then p=—q follows 
y 2 2 
and lo. = = = “ F 2) It seems that a theory in which a uni- 
form s.d. is impossible cannot be considered as a general linkage 
theory for tetraploids. In this theory, which has been designed in 
order to explain a particular set of data from observations, the prob- 
lem of tetraploid linkage is reduced to diploid linkage. 

If the theory is generalized to s=2 , m=8 it can be seen that of 
the 216 possible gametes, certain types do not appear. Hence this 
theory may apply to a situation in which certain classes of gametes 
are absent. . 


4. The General Segregation Distribution. 


We now consider a 2s-ploid with m loci each determined by one 
of r alleles. In the formation of a new individual corresponding to 
each locus, a set of s values is transmitted. We shall define a gen- 
eral chromosome segregation distribution. Denote by 


( i VV (C2 aeons (i; SoA ot), 


the probability that the transmitted set of s values which corresponds 
to the ith locus contains «, maternal values (and, consequently, s — «; 
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paternal ones) where ¢; = 0, 1,---- 8;¢1=1,2,----m. Each single 
A(é: €2 ---+ &m) is the probability for a specified set of maternal and 
of paternal values and 


(2) 


is the probability for a specified set of maternal values (see section 
2). We then have 


On..s 0....8 8 2 8 2 
ps ( )| | )| A(é. +++ Em) 
£4 Em €y Em 


| : ) ey ( B ) Mere) =] 


1 Em 


ie ) ey) at (eee) (28) 


S 


(29) 


As in section 2 of this paper and in earlier papers (Geiringer, 1944; 
1948a) we have to assume a symmetry relation 


ieere. yd (Sera 8 B,) (30) 


We see that for s=1 and general ““m” (27) reduces to a linkage 
distribution (l.d.) as defined and used in H. Geiringer (1944), and 
for m=1 and general s to the s.d. defined and used in H. Geiringer 
(1948a). The s.d. (27) is a discontinuous (arithmetical) m-variate 
distribution, each variable equal to 0, 1, ---- or s, and for which 
the relations (29) and (30) hold. 

We next introduce the marginal distributions of this s.d. We 
shall show that by an appropriate “marginal summation” we are able 
to reduce the “‘s”, the “m’, or both. From I(e, ---- &m) we derive by a 


(vy, +--+» + y,)-tuple summation, (0 = vi S 8), the marginal dis- 
tribution, 
Serato en AE? (€, Pade o Bai) 
(31) 
O..0:V4 Ose Vm v1 ) ( =) 
am a ae l cea EE Gee ee Sar ) 
2 "3 & ; 91 1m 
where | 
ee=0, 1,2 (8-11), G1, m). (32) 


In a symbolical way which is both useful and suggestive, the defini- 
tion (31) may be written 
[0¥3).--- (Ym) (gp, +e. Em) ~ 1) (e,) «++: LY) (en), (33) 


where J (e,) is the marginal distribution introduced in H. Geiringer 
(1948a, equation (32)), with », and «, standing for vy and a used there. 
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In fact for m1 and with » and « instead of v, and «,, the definition 
(31) gives 


1) (a) =| ‘ )u@) +( ' Jua+1) to +( ‘ )ua+»). (34) 


For s=m=2 the definitions reduce to those of section 2 of the pres- 
ent paper with the marginal distributions defined in (10), (11), and 
(13). For the distribution (31) we have 


0....S—V4 0....8—Vm sS—y S— Yn 
Sue Ss [meta (6, 00+ em) +f 4 Ly ( eae (35) 
€y Em 


&4 Em 


a relation that corresponds to equation (34) in H. Geiringer (1948a). 
The symmetry relations for the marginal distribution (31) follow 
from (30): 


[¥1) +. Vm) (e; eds Em) = [ (11)... (Vm) (s a SH pe, SD, Em) , (36) 


For », = +--+» = ym, = 8s —1, (86) reduces to the symmetry rela- 
tion (8’) of Geiringer (1944) and for m=—1 and general “s” to the 
symmetry relations (33) of Geiringer (1948a) (see also (29)). With 
these definitions the marginal distribution which corresponds to the 
normal organism, s==1 , becomes 


AS tae SOS Sa Brora a = hyo... (Er +++ Eis (e,=0 ’ 1c (37) 


We call it linkage distribution, while for s > 1 we speak of the segre- 
gation distribution. 

Another important particular case is obtained if we put », = 7. 
2+) = my = 0, vm = 8 in (31). This amounts to reducing the num- 
ber “m” to “m—1” without changing the “s” for the first (m—1) 
factors. The result is a distribution which no longer depends on’e;,? 


J(0)++-(0) (8) (g, ++ &m-q 0) 


0....8 8 rs (aa eat ie 
— DS Ue, -+ ++ Ema Feo) i ) ? ( ) ’ 
Nm 


Nm i=1,----, (m—1) 


(38) 


since, according to (32), em = 0. Finally, applying this type of sum- 
mation to (m—1) of the m characters we obtain the marginal dis- 
tributions of order one, and general “s” such as 


HO) (#)--01(8) (3, 109-90) = 1, (8,) a ayn (39) 


This is the s.d. for the first locus which has, because of (36), the 


required symmetry properties. There is of course such a distribu- 
tion for each character, such as 
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oO OO 4 on 0. 27.0 +--+, 0) 
40 
=I; (8;), (e,=0,----,8), S 
and 
L; (e;) =1;(s—s;). (41) 


We shall need all these marginal distributions in the proofs of our 
limit theorems. 

It is not difficult to write down any such marginal distribution 
if we make use of the “symbolical’’ notation introduced in (33). As 
an example consider for s=m=3: 


123) (000) ~ L® (0)1™ (0) 1 (0) 
== [(,:+,26 + lL]. [+4] - [4 +324 3b ':4+ &] 
= 1(000) + 31(001) + 31(002) + 7(003) 
+ 1(010) + 31(011) + 31(012) + 1(013) 
+ 21(100) + 61(101) + 61(102) + 27(103) 


(210) + 31(211) +°32(212) +1(218), 
or, in the same way, 
123.076) (100) — fi, +22. +1,]- [+4] -[L 4+ 3.+ 3b 4+ 4]. 


We finally wish to define independence for our general s.d. We 
shall call [(e, ¢,---- ém), (e; =0,1,---- 8) independent if 
L(e,-+++ &n) = L (ex) 22+ bn (Em) - (42) 
Here, as we easily see, J;(e;) are the marginal distributions defined 
in (40). It then follows easily that if I(e. ---+ &m) resolves mto a 
product, according to (42) the same is true for all marginal dis- 
tributions 


1Y) +) (ge, +--+ En), (e,=0 : il yore (s—yv)), 


and, eventually, 


1 
| (e, st Em) = baton (é€, yee Em) = Qn" (43) 
for all combinations (e,, +--+ &m) where s, = 0,1. It is well known 


that (43) represents G. Mendel’s original assumption of “indepen- 
dent assortment.” For m=1, general “s’”, the equation correspond- 
ing to (43) is equation (38) in H. Geiringer -(1948a). Obviously 


(42) is satisfied if the ( 2 yy values of the s.d. are all equal to 
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25a : 
each other; hence each equal to 1 y ( - )] . In this case each 


Lig. ' (8, -** Em), (8s = 0, 1) of the Ld: equals 1/2” and each 1; («i), 
2s 
Ss 


(e;,— Otel, 1S) equals 1 /( IE 


5. Recurrence Formulas for m=s=2. 


It seems useful to investigate the recurrence relations for the 
particular case m=s=2 by explicit computation before dealing with 
the general case in which explicit computation is not possible. We 
shall denote, as before, by p™ (A*a?*B“b?) the probability dis- 
tribution of gametes in the nth generation with sum equal to one ac- 
cording to (1). From this distribution we derive two marginal dis- 
tributions of order three. Using the variables oa; equal to A or toa, 
and £6; equal to B or tob, (t= 1, 2), we introduce 


D123"”) (a1026,) = =e ape (1,021 Ps) = Py24" (a:02he) . 
Bo 


We may omit the subscripts and have the two distributions 


p (a,a26,), and p” (afife), (44) 
with 
Son (a1); ) =a Loe pp (a,6:B2) =1. 
Q@ a2 Bx a1 Bi Be 
There are three marginal distributions of order two, 
po (a,a2) ’ pm (a,8,) ’ pm (B,B2) ’ (45) 
where, e.g., 
p™ (ai02) = 5 p™ (a,026;) =TDTV™ (arOe08,f.) ; 
Ba Bi Be 
with | 


22 (a2) =1, Sd p™(a6:) 15 ete; 


a, a2 ay Bx 


or more explicitly, for example, 
p (A?) =p (A°B) + p (A) 


=p (A2B?) + 2p™ (A?Bb) + p™ (A2b?2) ; 
and, finally, two marginal distributions of order one, 


p™ (a), pe (6), (46) 
where 
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DING) —= Sip (o.d,) =>) 9™.(0,8,) S425 
a2 Bi 
a a as 3 Do (05028382) , 


a2 Bi Be 
and 


Dp (az) = 1; Sp (By) == 15 
ay Bx 


We now wish to derive p(a:02$:62) or, more specifically, 
p‘”) (AB?) from the distribution of genotypes, w™, in the nth gen- 
eration. For this purpose we shall consider each one of the various 
terms in (2) and multiply it by the corresponding probability values 
of the s.d. For the sake of brevity let us now write p™ = p, 
w™ =w,p™ =p’, w™) = w’. Consider, for example, in accord- 
ance with (2) and (8), the term 8w(A?B?; AaBb). If we derive 
p(A?B?) this expression is multiplied by (Ass + 2dy + 2do, + 4d41) 
where we write Jag for A4(a, 8). In fact, doo is the probability that 
the four values AABB in p(A?B?) are all maternal values; 2/,2 is the 
probability that of the two A-values one is a maternal, the other a 
paternal one; while both B-values are maternal heritage, and simi- 
larly for 22.,; finally, 44,, is the probability that one A-value will be 
paternal, the other maternal and the same holds for the B-values. On 
the other hand, for example, the coefficient of w(a?b?; A?Bb) will be 
zero since (under chromosome segregation) the four values AABB 
cannot be transmitted by an individual of this type. If we consider 
each w-probability in this way we find 


p (A?B?) 
= w(A2B?; A2B?) - 1+ 2w(A2B?; A?b?) (Aes + die +Aoz) 
+ 2w(A2B?; a?B?) - (Aso + 4dor + doo) + 2w(A?B?; @?b?) doo 
+ 4w(A2B?; A?Bb) - (Aes + 4dre + doe + 2dor i+ 8Air + 2hor) 
: (47) 
ii 
+ 8w(A2B?; AaBb) - (Ase + 2dre + 2der + 4dr) Teo: 
+ 16w(A?2Bb; AaBb) - (ao + 2d) +°:°: 
+ 16w(AQBb; AaBb) -i4,.. 


There are on the right-hand side of (47) twenty w-values of which 
eight are written. Note that the coefficients of each w-value are noth- 
ing but marginal probabilities of our s.d. (see section 2). For ex- 
ample, Joo + 4dw + dco = (2), the coefficient of 8w(A?B?;AaBb) 
which has been discussed in detail, (Aiz2 + 242 + 241 + 41), equals 
l,.(11), while the coefficient of 4w(A?B*;A*Bb) is (es + Bie + Ace 
Re ee Oi DA et (11) be (01) = 1/2. It would not be 
difficult to describe the general structure of this type of formula. 
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Now we write the same formula for p’ and w’. Next, we express 
w' = w in terms of the » = p\ by means of (4). We obtain a 
formula similar to (47), 


p' (A2B*) = p (A®B*) p (A?B*) + 2p (AB?) p(A%*) 
- (Ags se AW + Ries) a oes eyes te 


The decisive step is then to rearrange this formula by collecting 
on the right-hand side those p-products which are multiplied by the 
same jag. Before doing this the coefficient “1” of p(A2B?) p(A?B?) 
is to be replaced by the sum of all A-values, using (9). It then turns 
out that the coefficient of each dag is the product of two marginal 
probabilities as introduced in (44)—(46). (Note that while before 
the coefficient of each w-value was a marginal probability of the 
s.d., now the coefficient of each 4-value is the product of two marginal 
probabilities of the p-distribution.) 

The final result is the following simple formula where the d’s 
have been replaced by the respective -values in consideration of (5) 
and (8), 


p' (A*B*) = p(A*B?) [1(00) + 1(22)] 
+ 2p(A) p(AB?) [1(10) + 1(12)] 
+ 2p(B) p(A2B) [1(01) + (21)] (48) 
+ 4p(AB) p(AB) -1(11) 
+ p (A?) p(B) - [1(20) + 1(02)]. 


This formula can be interpreted in the following way: 


woo-(2) (3) 


is the probability that a maternal A-values and / maternal B-values 
will be transmitted as well as (s — a) paternal A-values and (s — §) 
paternal B-values; p(A‘°B8) is the probability of a maternal gamete 
With « A-values and § B-values; while p(A*-* B*-*) is the probability 
of a paternal gamete with (s—a) and (s—f) A- and B-values respec- 
tively, in the mth generation. Under random breeding the sum over 
all a and § of products of these three probabilities gives the prob- 
ability of a gamete of type (A?B’) in the (v+1)st generation. The 
corresponding formula for s=1, m=2 has been given by H. S. Jen- 
nings (1917) and the formula for s=2, m=1 by J. B. S. Haldane 
(1930). If for s > 1, we consider chromatid segregation, the recur- 
rence formulas for m=—1 and for m=2 no longer show this particu- 
larly simple structure. 


In an analogous way a formula for p’(AaBb) may be derived. 
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It is easily seen that the structure of these formulas does not 
depend on the number of alleles. If there are more than two alleles 
for each factor, there will be more terms in the formula which takes 
the place of (47). However, in the decisive transformation where 
we determine the coefficient of each dig, the sum of all p-products 
which multiply the same dag will, exactly as before, resolve into the 
product of the same two marginal probabilities as before since each 
marginal probability is then defined by summations over r rather 
than over two terms. The formula (48) and the following formula 
(49) thus hold with the understanding that A , a as well as B , and b 
are just any two out of an arbitrary number, 7, of alleles. The re- 
sult is: 


p' (AaBb) = p(AaBb) [1(00) +1(22)] 
+ [p(A) p(@Bb) + p(a)p(ABb)][1(10) + 1(12)] 
+ [p(B) p(Aab) + p(b) p(AaB)][1(01) +1(21)] (49) 
+ [p(AB) p(ab) + p(Ab) p(aB)] - 21(11) 
+ p(Aa) p(Bb) - [1(20) + 1(02)]. 


From (49) we obtain (48) by putting A=a, B=b. 

Thus the distribution of gametes can easily be computed from 
generation to generation. The distribution of zygotes then follows 
for any desired generation by means of (4). 

Obviously these formulas are valid for any particular choice of 
the parameters entering in the definition of the s.d. If for example 
we consider random chromosome segregation (each J-value equals 
1/36) we obtain: 


36p' = 2[p(AaBb) + p(Aa) p(Bb)] 
+ 4[p(A) p(aBb) + p(a) p(ABb) 
+ p(B) p(Aab) + p(b) p(AaB)] 
+ 8[p(AB) p(ab) + p(Ab) p(aB)], 


or, for Haldane’s theory, the five coefficients in (49) are to be re- 
placed in turn by 92/3, pq/3, pq/3, (p? + @ +pq) /3 , and p?/3. 

In the second part of this paper we shall establish an analogous 
recurrence formula for any “m” and “s”. By means of these recur- 
rence formulas we shall investigate the equilibrium problem and the 
limit behavior for our basic distributions, the distribution of gametes 
and the distribution of zygotes. 
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Steady state kinetic relations previously developed for closed, 
homogeneous systems are extended to open systems consisting of geo- 
metrically confined regions of arbitrary shape. The generalized system 
of consecutive reactions is considered to occur within the cell, and the 
cell plus environment are treated as an open system. The diffusion con- 
dition is imposed upon the kinetic solution for various special cases and 
the method of general solution when all products, reactants and inter- 
mediates, may enter and leave the cell is indicated. 


In a previous paper (Hearon, 1949), hereinafter referred to as 
I, the system of consecutive reactions 


Ce Oy Ag Dy 
> 

Ce AS SX, PB 
=> 


(1) 


C; ie XG << xi in B; 
— 


Ca Bis ila ey B n 
was considered to be in the steady state. As noted in I, this assump- 
tion requires that te = 0 for all 7. The vanishing of the time 


derivative of all [X;] in turn requires that the concentrations of all 
reactants [C;],i1—=0,1,---- n and all products [Bi], ee a Ee 
---. 2 must be time independent. (Exceptionally, this restriction may 
be removed from [B,] if k., = 0.) These conditions may be achieved 
in several ways. All the ways discussed in I applied strictly to the 
“closed systems” considered there. When in vivo systems are consid- 
ered in biology, it is most usual to deal with open systems. In such 
cases the reactants and products may be maintained constant by sup- 


+This work was done while the author was a Senior Research Fellow of the 
National Institute of Health. 


83 


84 KINETICS OF BIOLOGICAL SYSTEMS 


ply to the system from the environment and flow from the system. 
The situation is then that the system proper is the set of reactions 
being considered; the reactants are maintained constant by sources 
in the environment and the products flow to sinks in the environment. 
These sources and sinks must now have the property that addition or 
withdrawal of material does not change their levels. This require- 
ment is met physically by having the external concentrations con- 
stant by adherence to short experimental times, or by employment 
of environment of such size that quite regardless of the absolute 
levels (concentrations) the addition or withdrawal of even an ap- 
preciable amount of material changes the levels only negligibly. 

Of particular interest is the flow of material from the ambient 
medium into a cell, its conversion there through a set of reactions 
to products and the flow of products into the medium. In general the 
reactions will occur only within the cellf and the cell plus environ- 
ment constitutes an open system into the description of which must 
enter the parameters governing the flux of matter into and out of the 
cell. The rate of consumption of any C; will depend upon the con- 
centration of C; in the system and upon the concentration of every 
other component in the system. The concentration of C; in the sys- 
tem in turn depends upon the rate of consumption according to the 
equation (Rashevsky, 1948) 


[Ci] = [Cio — Ao, Vo,» (2) 


where [Ci], is the external concentration of C;, Ve, is the rate of 
consumption of C; and Ae, is, as defined by N. Rashevsky (1940), 
the total diffusion resistance of the cell to C;. In general Ag is a 


function of the diffusion coefficients of C; in the medium and in the 
cell, of the permeability of the cell membrane to C; and of the major 
and minor axes of the cell (Rashevsky, 1948). Exactly similar con- 
siderations apply to any product B; and the equation analogous to 
(2) is 

[Bi] = [Bilo + As, 0s, (3) 


While many systems of this type have been solved (see Rashevsky, 
1948) and indeed some rather complicated cases (Landahl, 1939), it 
seems desirable to attempt to lay down some general formulations 
based upon the steady state solutions previously discussed (I). The di- 
rect procedure, which has been previously used, is to write v; , equiva- 
lent to the net rate of an individual step in the system (1), for every 


This statement already implies that certain component i 
involved, may not leave the eal” sense mise set 
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C,, X; and B; and to substitute into these relations the correspond- 
ing equations (2) and (3) and solve the set for the desired rate. 
There will be considerable advantage in making use of the steady 
state solution and solving this with the diffusion equations (2) and 
(3), since the assumption is always made that the reactions in the 
cell are in the steady state. This procedure, when it can be used, has 
the feature that the expression for a given rate will be obtained in 
a form analogous to the expression for the homogeneous system, thus 
the effects of diffusion, geometry and size are more readily assessed. 
This will be illustrated first by solving the simplest possible enzyme 
system. 


The Simple Enzyme System. 
Let the mechanism be 
Bee Sh S > h+ op; 
—= 
using the notation of I, and let S and p only leave and enter the cell. 
The steady state solution is then (equation (10), I) 
_ ke[Eo][S]_ Va[S] 


a ea? (4) 
Kas] Ko 1 | 
where 
* _katkh 
a 
Equation (2) becomes 
[S] = [S]o— Av, (5) 
and substitution of (4) into (5) and solving for [S] gives 
S. One 
[S] ee 1/2V ([S]o —@)? + 4Ks[S]o, (6) 
wherea—=K + AV,,. Using (6) in (4) gives 
Val{Slo— Av) 
qd os 
Kg + ([S] os 


a eae: + 1/2/(STo— 4)? + IKASTo| 


ee a SS nn ln nnn 
A a _ eet 
Ks + joes + 1/2V (Sle — 4)? + 1KS1;| 
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Here and elsewhere we will denote by q the diffusion conditioned 
rate and by v the simple kinetic solution such as (4). From (7) it 
is seen that if [S]. is large, 


VnLS]o 
Ks + [S]o 


which is identical with v from (4) with [S] = [S].. From (4) and 
(8) it follows that the limiting velocity V», is the same for g and v. 
The conditions for (8) may also be that A is small (D; large, h large) 
or that v is inherently small (k. small or [E], small). Equation (8) 
then becomes the analytical statement of the fact that if the total 
permeability, 4 = 1/A (Rashevsky, 1940), is large or if the rate of 
consumption is small, then [S] ~ [S]. and v and q do not differ. 

Conversely, if A is large or if v is inherently large, then from 
(5), [S] ~ 0, and 


q= (8) 


q=Al[S]o- (9) 


In this case diffusion is completely limiting in the sense that the 
observed q is simply the rate of diffusion into the cell and does not 
contain the parameters, K and V,,, of the system proper. 

An intermediate case which exhibits a point which should be 
noted is: If [S] is small (K >> [S]), then from (4) and (5) 


TEE 
[S] Se eat pres (10) 
1+— 
K 
and (10) into (4) gives 
Vn Halo 
Pua, ( AVin ) Gh 
1+ — 
K 


From (11) it is clear that if either A or Vmn(=k.[Eo]) becomes large, 
(11) reduces to (9), but as A continues to increase g — 0 while the 
limit of gas V,, > o is given by (9). The above discussion is given 
in some detail because analogous cases will arise in the more general 
systems to be discussed. Finally it will be indicated how the para- 
meters K and A may be determined in this simple case. From (4) it 
is seen that if [S] = K, then v = 3 V,, and that the initial slope of 
the v vs. [S] curve is given by 


( aa ) [sI<0 ==. ve 
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From (5) when v = 4 Vn the required external concentration, 
[S].°/), is given by 


AVin 
[S$] 0/2 = K + a (13) 


Thus while v and q exhibit the same limiting velocity V,,, v reaches 
+ Vm at [S] = K and q reaches + V,, at 


AV on 
[S],»=K as 2 ? 


ie., the external concentration [S]) must exceed the internal concen- 
‘AV: 
tration [S] by the amount Ey The initial slope of the q vs. [S]o 


curve is 


é Vin 
( ) So eae (14) 
OLS] [8] =0 K+ AVin 


From (13) and (14) the values of K and A are obtained as 
cei 2 [S]o°/? mp ae Vin. 
= a : 
1 S] 2 
a=2( eh 


Mo Vg 


K 
(15) 


The above analysis supposes that V,, may be observed directly, 
i.e., when gq becomes independent of [S]. at large values of [S]» it 
is assumed that this value of g is Vm. In practice this assumption 
is not always a good one even for isolated enzyme systems (see, for 
example, Lipmann, 1939). 

The General System with C, Diffusible. 

In this case consider that only C, diffuses into the cell from out- 
side; all other components of the system are confined to the cell ex- 
cept the final product B,. If the final step is irreversible, ken = 0, 
it is immaterial that B,, is diffusible. The condition (2) which is here 


[C,] = [Ci]o— Ac,U (16) 
is easily included in the problem by writing the first member of the 
set (2) from I as 


[hose Kho, [Go] = Wy) — W101 » (17) 
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where w,° = k,[C.] [Ci]. and the a’s are as defined in I. Equation 
(17) results from the substitution of (16) into the definition of wi 
in (1) of I. The problem may also be solved more directly by writing 
the general solution of the system (see equation (3) of I) as 


w, IT wi 
1=2 
oh eee (18) 
> 3 
where db = WeWs **** Wy + WaWs-00* Wr He Ss = DL Wag es Weep is 


the denominator of (3) of I and w, = k,[Ci] [Co] . From (16) and 
(18) 
[Clo 


[Gal eee 19 
1+ £/¢ oe 
where 
B — Ag Ka [Co] W QZ. 
Substituting (19) into (18) gives . 
Jes[Ca] [Calo IE wi 
g a V2 
(1 + B/$)¢ 
(20) 


ee WWW; ot! Wy 
(1+h Ac [Co] ) We Wer s+ Wat W Wye Wyte PW 2 Won) 5 


Equation (20) differs from equation (18) only in that the first term 
in ¢ has been multiplied by 1 + ky Ac [Co] . It will be noted from (18) 


and (20) that while v is linear in [Co] , g exhibits a limiting velocity 
when [C,] is increased indefinitely. This limiting velocity is Ac [Ci]o- 


The interpretation is clear and is in accord with the general prin- 
ciples set forth in I. When [C,] is increased, the rate of the first 
step, and hence of the entire system, cannot exceed the rate at which 
C, can diffuse into the cell; the rate is Jo [Cr ]o. The dependence of 


q on C;,,7i > 2, differs in an obvious manner from that of v. For 
example, the limit of q when [C.] > o is 


w,? 
Let kyAc [Co] ; 


while the limit of v is simply w,. 
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The General System with C, Diffusible; i~0O,1. 


The same considerations apply here as in the above case, viz., 
only a particular C; diffuses into the cell and the final product B, 
diffuses out; it is again assumed that k_,,=0. The general solution 
for v from (18) may be written as 


naa , (21) 


where the T; are the additive terms of ¢; there will be n such terms, 
each the product of n — 1 w’s. The terms 7, through T;., will con- 
tain w; as a factor. Thus (21) may be written as 


n i t=-1 Tz 
IT wi >» 
Gee 2 = W; 
yee ee (22) 
n %-1 € 
Ww; + ST; =— 
= Ey Wj 
Since the terms T, through 7;_, will always contain WisWiis +++ Wn aS 


a factor, it can be shown that the denominator of (22) is 


WW 2+ Win 


i : (23) 
(WeWe Wi) (WW * + Wi) He Ft WaW_2W_o + + -W_ (5-2) 

The expression in the brackets is just Vi,:i as defined in I. For rea- 

sons that will become obvious it will be denoted here by V»,“. 

Then 


Ve" [Ci] 
y= ———; 
[Ci] + @; 
(24) 
bagi 
es &i 
Bhi pews ay 
ki > Te/wi 
&1 


It is now obvious that V,,‘“ is the limiting velocity for C; and may 
be obtained at once from (21). The general expression ®; can be ob- 
tained by deleting [C;] from the sum of those P terms which contains 
it as a factor and dividing the sum of the remaining terms by k; times 
the result. With v in the form (24) the condition [Ci] = [Ci]o — Ac,Y 
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may be straightforwardly imposed. The solution of 


Ao Vn? [Ca] 


.1= [¢C.], -———— (25) 
[Ci] = [Ci]o iG +o, 


[Cilo ay (9%; se Ao Vn?) 


[el 9 
(26) 
V{{Ci] aes (; sis Nae Vim) fe rE A[Cilo ; 
PRAM Mn Ree EOL RE 
as a, 
Using (26) in (24) gives 

Vn) {4+ H 
26,+ {G+ H} 


where G = [Ci]lo — (®& + Ao Vn) and H = VG + 4[Ci].®:.- 
This form may be briefly discussed as follows: When [Gale 0s 


g > Vn'® as expected from the mode of derivation of (27). When 
ho, is small, gq becomes 


SVG OLE Ne 
oe [Ci], + %; : 


which is formally identical with the same situation for the simple 
case (equation (8)). When ho, becomes large, 


(28) 


q=4e, [Clos © (29) 


identical with (9) for the simple case. As in the simple case there 
is an intermediate situation analogous to (11). 


When the dependence of q upon [C;] is considered, one must 
distinguish between the cases in which k >i andk <7. 


When [C;] , k > 7, is increased indefinitely, the only quantity in 
q which varies is ®;. It is easily shown that ®; approaches the 4; 
for the first k — 1 steps of the system and q approaches the corre- 
sponding q. Since 
ke 1 Tio Ti aT ets 
kw; i a 4A ss Se 8. ate Aa 


and only the terms T, through 7;., contain [C;] , 


i 


(30) 
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tee Egeland 
Lim #; = a (31) 
iglaa kw; ra T, +++: Ti4 ; 
Both numerator and denominator of (31) contain ww +--+: W_ as 


a factor, and comparison of (30) and (31) shows that the limit (31) 
is Just 6; for a system of k — 1 steps. If ,™ is ®; as defined for a 
system of 1 steps, then 


Lim ©, = 6;*1) , cee (32) 

[Cz] 00 
When [Cx] , k <i, k #1 is increased indefinitely, 6; becomes 
zero as is clear from (30), and from (23) V,,‘°) becomes V»(°). These 
particular values in (27) give the corresponding g. If k= 1 and [Ci] 
becomes indefinitely large, from (27) q takes the indeterminate form 


ao -0 : 

26,;™ : 
for from (26) [C;] becomes zero and from (23) Vn? = o. It is 
easily seen, however, that from [C;] = [C;], — Ao,v the actual limit 
is 

q = de [Ci]o- (33) 
In summary then 
Ling = V6? 


[Cilo>oo 


CG VP) ak oe 
Lim q(@i, Pn) | g(@™=0, Vn) k<i (84) 
ie k#1 

det Cals PS ihe 


The above cases, C, diffusible, C;, 1 # 1 diffusible, and the simple 
enzyme system are of somewhat restricted applicability in that many 
common systems are such that several C; and several B; (produced 
in a reversible step) are diffusible. It should be noted, however, that 
if a single diffusible substrate is being utilized through a typical dehy- 
drogenase system as discussed in I, this corresponds to the case C,; 
diffusible; for the enzyme, the complex E - S and the high molecular 
weight carrier, (such as flavoprotein) are certainly confined to the 
cell. The only requirements are that the permeability to the coenzyme 
be such that its flow from the cell is negligible and that the product 
is produced in an irreversible step. In practice these are frequently 
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good assumptions. Similarly if certain reactions precede the step 
in which the diffusible substrate enters, this corresponds to the case 
of any C; diffusible. The case of enzyme systems with dissociable 
prosthetic groups treated in I may be handled as the simple enzyme 
system, for equation (46) of I may be put in a form which is formally 
identical with equation (4) of this paper. The method of solution is 
identical with that given here and it is only necessary to replace Vin 
and K in equations (4), (6), (7), (10) and (11) by the appropriate 
quantities gotten from (46) of I. 


Any C, and any By, Diffusible. 


In the general solution for v , [Ci] and [Bx] occur only in w; and 
w_, and v may be put in a convenient form for simultaneous solution 
with equations (2) and (3). The form depends so strongly upon the 
position in the system of C; and B; that only the method of solution 
will be indicated. Consider a system of four steps with k..=0. The 
expression for v is 


W1W2W3W4 
= : (35) 
WoW3W, + W1W3W4 + WaW Ws + W1W.W_s 
If only C, and B, are diffusible then (35) may be written 
[Ci] L 
v= : (36) 
M + [B.]N 
Substitution of (37) 
ral a [C.] Otme 1 Ac V 
(37) 


[B.] = [B2}o + Apu 


into (36) gives a quadratic in v which may be solved for v in terms 
of the A’s and external concentrations of C, and B,. It is clear that 
many combinations of B’s and C’s diffusible may be solved in this: 
manner. The manner in which the positions of B; and C; in the sys- - 
tem influence the degree of the equation to be solved raises an inter- . 
esting question. The answer in a given case is immediately obtained 
from inspection of (35); as yet no general answer has been obtained. . 
For example, the following table gives the degree of the equation in 
» for all possible combinations in the system under discussion when 
any C; and B,, are diffusible. 
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a k degree 

it Te 2.05 quadratic 
2 ” 9 

3 1 cubic 

3 aOLE quadratic 
4 Or 2 cubic 

4 3 quadratic 


Similar considerations apply when more than one C and one or more 
B are diffusible. 


The General System, X, Diffusible. 


Again the perfectly general case has not been worked out but 
the method will be clear from the example given below. 
If X, is diffusible consider the steady state requirement 


a[X,] | 
dt = W, — W4[Xi] — Wie [Xi] + we, X.—T(X1) =0, (88) 
where I'(X,) is the flow of X, from the cell. From (38) 
W, — W(X] = wo[X.] — w.[X.] + (NG). (39) 


Also in the steady state 
VO Wy Wa XG | 
= W,|X,| — wel Xo] +F(X,) 
=w,[X.] — w-s[Xs] (40) 
Ny [i Wy 8 
The quantity I’(X,) is now given by 
I'(X,) See (Xa —= [Xx16); (41) 


where ie ies aie and [X,]. is the external concentration of X,. Tak- 
ing the case in which n = 3 and k., =0, (40) becomes 

Di 0g ay 
1L= (ws + dx )ay— Wat — Ax, [XiJoot ce 
1= W530. , 


which is the set corresponding to (2) in I; the notation used here 
is the same as that in I. The solution of (42) for v = 1/a which in 
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our notation will be called q is 
W (We a Ax,)Ws — WAW3 [X,] oAx, 


qd = ‘ (43) 
(We + Ay) Ws + WsWs + WaW-2 


Equation (43) differs from the corresponding expression for v (ob- 
tained from (48) by setting Az = 0) in that w, is replaced by w. + Ax, 


and the subtractive term in [X1]. appears in the numerator. Simi- 
larly, if C, and X, are diffusible the set corresponding to (42) becomes 


1 ae Kyo | [Co] — WEA ae W_104 
i Woy — W2Ae (44) 
1 (We Ax. )o5— Ax, [Xs Joa, 
where the first member of the set was obtained by using equation 
(17). The solution of (44) is 


WWe (Ws + Ax.) a W-1W-ohx, [Xe]o 


q (45) 


(1+ Ky[Co] Ac.) (ws + dx,) We + Wa (Ws + dx,) + Wawa 


It is readily seen that if all of the X’s are diffusible, the set (40) 
becomes 
Vv=W,— wWa[X,] 
= w,[X,] — w+ PAY ay T'(X,) 
= Ws[Xo] — ws[Xs] + P (Xe) + I'(X3) (46) 


=WrlXral — Wn + 3r(X) a 


4=1 


It is also obvious that if in addition to any number of X’s one 
or more C’s and B’s are considered diffusible, the equations (2) and 
(3) may be substituted into the corresponding w’s in (46). The so- 
lution of (46) then gives the perfectly general case. From (46) it is 
evident that only the first member will be linear under all conditions 
but the set can in principle be solved systematically; the exact meth- 
od will depend upon the case at hand. It should be noted that the dif- 
fusibility of the X’s does not influence the linearity of the members 
of (46); the diffusibility of any B and C except C, causes the non- 
linearity of the corresponding member of the set (46). 
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: When two polished metal spheres in contact with each other are 

» immersed in a corrosive solution, their surfaces gradually corrode, leav- 
ing, however, a non-corroded zone around the point of contact. The size 
of this zone is a function of the size of the spheres, concentration of the 
corrosive, and of the time. The phenomenon has been first observed by 
F. C. Besic on human teeth immersed in hydrochloric acid, and studied 
on metal balls. In the present paper it is shown that such a phenome- 
non may be due to gradients of concentration of the corrosive in the 
neighborhood of the point of contact, caused by the chemical reaction 
which consumes the corrosive. Approximate expressions for the size of 
the uncorroded zones as a function of the size of the balls and as a 
function of time are derived and found to be in fair agreement with 
F. C. Besic’s data. 


It has been observed by F. C. Besic (unpublished) that when 
extracted teeth were immersed in various acid solutions, all of them 
showed an area one to two millimeters in diameter that was devoid of 
corrosion around the point where the crown of the tooth touched the 
beaker in which it was placed. A similar phenomenon was noted 
when two teeth were placed in contact with each other and immersed 
in acid. In all cases the uncorroded area was definitely larger than the 
area of possible contact. 

It has been suggested that the phenomenon was due to a sufficient 
drop of concentration of the acid in the narrow space between the 
teeth, around the point of contact. The acid was used up by the corro- 
sion reaction on the surface of the teeth. Because of the narrowness 
of the space in the immediate neighborhood of the point of contact, 
the amount of acid diffusing into that space may be insufficient to 
cover the shortage created by the reaction. 

To investigate the phenomenon under better controlled geometri- 
cal and physicochemical conditions, the use of steel balls of different 
diameters in contact with each other was suggested. Hydrochloric 
acid of different concentrations was used as a corrosive. In this way 
F. C. Besic investigated the variation of the size of the uncorroded 
zone as a function of the concentration of the acid, the time, and the 


size of the balls. 
A separate study of the rate of corrosion of steel by hydrochloric 
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acid was also made by F. C. Besic for the purpose of the present 
study. The result indicates that the rate, q, of consumption of HCl 
expressed in gm cm? sec”, when plotted against the concentration, 
c, of the acid, gives a curve convex downward. The accuracy of the 
measurements was insufficient to establish exactly the form of the 
function 


q=f(c). (1) 

In the present study we shall put, as a very rough approxima- 
tion, 

f(c) = ke’. (2) 


id 


Since the acid is consumed only at the surface of the steel balls, 
everywhere in the solution we have 


V7e=0. (3) 


Denoting by D the diffusion coefficient, and by » the normal to 
the surface, positive when directed outward, we have at the surface 
of the balls 


0c 
D—+ke’?=0. (4) 
Ov 


Considering the size of the beaker as very large compared to the size 
of the balls, a condition satisfied in the experiments of F. C. Besic, we 
have at infinity 

¢ = c¢, = (const). (5) 


Because of the nonlinearity of the boundary condition (4), and 
because of the peculiar shape of our boundary, an exact solution of 
the system (3)—(5) is not to be hoped for. The use of the approxi- 
mation method (Rashevsky, 1948) is indicated. 

Let the two balls be of the same radius, R. Using the notations 
of Figure 1, we have 
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ere es i x=Ksin ¢; sng=x2/Rk. (6) 
y= R(1—cos¢). (7) 
The area, S , of the cylinder of radius x and heighth 2y is equal, be- 
cause of (6) and (7), to 
S = 4axy = 4nR? sin 6 (1 — cos 4). (8) 
Equations (6) and (7) also give 


== —_ =e 
y (1 a2 =). (9) 


If x/F is small, as is actually the case in the uncorroded region, we 
may expand the radical, keeping only the lowest power terms. This 
gives 


 icstlee ae (10) 


Approximation (1) holds fairly well up to x/R = i, since #!/R* = 
1/16. As a very rough approximation, we shall assume it to hold in 
the whole range,O <a <R. 

The element ds of area of the spherical zone between ¢ and d¢ 
is, because of equation (6): 


2naxdx 


CS — oasis oo ——— (11) 
V1l—27/F 
For sufficiently small values of «/R, we have approximately 
ds=2nxdz, (12) 
or, 
S=nx. (13) 


Equations (8) and (10) give 
S27 kK . (14) 


Consider the space limited by the two spheres and the cylinder 
(Figure 1, shaded). If é denotes the average concentration in that 
space, then, neglecting possible gradients in the direction of the axis, 
00’, we have for the total rate of consumption of the acid, from (2) 
and (13) 

Qn x? ke. (15) 


If c is the concentration at the surface of the cylinder (also con- 
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sidered as constant at all points of that surface), then the average 
gradient is 


2(e—6é)/a, (16) 
and the total inflow, because of (14), 
2DS (c= 6) /x = 4a Dx? (c— €)/R. (17) 


In astationary state (15) and (17) must be equal. Equating (15) 
and the right-hand side of (17), we find: 


kR@ + 2Dé—2De=0, (18) 


om 


(/1 + 2kRe/D — 1)D 


kk (9) 


C= 
The other root of (18) is negative, and, therefore, has no physical 
meaning. 

To obtain a second approximation, we shall now set the exact av- 
erage inflow at x, that is, DS(dc/dx) equal to the average consump- 
tion inside of the considered volume, 2a2?ké?. The first expression is 
equal, because of (14), to 


2nDx* de (20) 
R da 
The second expression is equal, because of (19), to 
2a? (\/ Ve Behe LD 1) 7 
Re ; (21) 
Equating (20) to (21) we find: 
de D(\V1+ 2kRe/D—1)? 5 
dx kRx : oO 
or, by denoting by K an integration constant: 
de D 
ee eee (23) 
(V1+ 2kRce/D —1)? kR . 


The integral in (23) can be readily evaluated by using the sub- 
stitution 
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After simple rearrangements, we finally obtain, with A = exp. 
KkR 


D 


1 


= A(V14 2kRe/D—1)e virmEo. , (25) 


When x = R, the concentration, c, is practically equal to ¢; 
therefore, A is evaluated by the requirement that forc=c,2=R. 
This gives finally: 


se al 1 
— Pp AYE IS eg © VEGRR/D-1 yis2kRe/D-1 (26) 


Apparently a minimum corrosion of steel is necessary in order 
to become noticeable. For a given time, t, we need a given concen- 
tration c* to achieve this. Since the rate of corrosion is approxi- 
mately proportional to c*?, approximately, 


c*2t = o? (= const). (27) 


Hence, the radius, 2, , of the uncorroded zone is obtained by put- 
ting in (26): 


poo /\/ i: (28) 
We thus obtain 


1 Z 


(it 2 Ra/Dijt— eae giana “exons 
R ei a ees @ VintkRe/D-1 VayzkRa/dyi-a + (29) 


1 + 2kRe/D — 1 
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For fixed values of c and t, % varies approximately linearly 
with R. This is borne out by experiments as can be seen from Figure 
2. The latter was obtained from graph 3 of F. C. Besic’s paper by 
plotting the values of the diameter of the zone against the values of 
R for three different relative values of the concentration. 

In alt of F. C. Besic’s experiments the uncorroded zone consisted 
of two concentric parts: the inner and the outer zone. The outer zone, 
while not appreciably corroded did not show the brilliant smoothness 
of the inner one. An explanation of this phenomenon apparently 
requires a more detailed analysis than the one given here. In line 
with the approximation assumed in (27), our %o refers to the radius 
of the outer zone. Empirically it is found that the radius of the two 
zones are approximately proportional to each other, and therefore 
the same relations hold roughly for both. As ¢ increases, 2o decreases. 
For very small values of t, 2% varies as 


at aS 
tore —— exp (—p/W1). (30) 
Wt 
For large values of t 
1 e 
By oc ——er!, (31) 
Vt 


The available data are not uniform enough to compare with the 
equation (29). They show, however, a general conformity with (31). 

For constant R and t, 2 decreases with ¢,, in qualitative agree- 
ment with experiments. For very large values of ¢ , 


il 


Ho & 


eve, (32) 


VC 


that is, x, tends to zero asymptotically. 
If, instead of (2), we assume a linear relation 


f(c) =pe, (33) 
then by the same method, we find 


to = RK ( : ee ie ~ (34) 
CoVt 
Putting 
PEt 2D 


2pR Us (35) 
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we have 


Lu 
log %) = const — oi log t+ log R. (36) 


xX (1) R = 1.225 
© (2) R =0.95 
oO (3) R =0.64 
A (4) R =0.32 


2X9 IN CM. 


—> t IN HOURS 
FIGURE 3 


Assuming that relation (33) holds approximately for small values of 
c (large pH), we may compare (36) with available data. Figure 3 
shows %p vs. t on semi-log paper for pH = 3, which is the next high- 
est used. This figure is drawn from data obtained from graph 5 of 
F. C. Besic’s paper, and refers to the inner zone since data for the 
outer zone are not available. In figure 3 we have R = 0.82, u= 04. 
Taking p = 3 X 10-* cm sec? and D = 2 X 10° cm’sec’, we find, 
from (35), «= 0.7, which is of the correct order of magnitude. 
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_Wken an individual grows up in a society, he learns certain be- 
havior patterns which are “accepted” by that society. He may in general 
have a tendency toward behavior patterns other than those which are 
“accepted”’ by the society. This tendency toward such unaccepted be- 
havior may be due to a process of cerebration which results in doubt 
as to the “correctness” of the accepted behavior. Thus, on the one hand, 
the individual learns to follow the accepted rules almost automatically; 
on the other hand, he may. tend to consciously break those rules. Using 
a neural circuit, suggested by H. D. Landahl in his theory of learning, 
a neurobiophysical interpretation of the above situation is outlined. 
Mathematical expressions are derived which describe the social behavior 
of an individual as a function of his age, social status, and some neuro- 
biophysical parameters. 


In our previous writings (for a summary see Rashevsky, 1948a) 
we pointed out that a logical extension of the mathematical. biology 
of behavior leads into mathematical sociology. The environmental 
parameters, which determine some of our reactions, are themselves 
determined by the reactions of other individuals. Thus we are led 
from the “problem of one individual” to the more general “problem 
of ” individuals,” and thence into the domain of social sciences. 


In our book (Rashevsky, 1948a) we attempted to construct a 
mathematical system of social sciences by starting with some plau- 
sible postulates about the interactions of two or more individuals. 
In principle those postulates should be derived from the equations of 
the mathematical biology of the central nervous system. As a tem- 
porary expedient, these postulates were introduced de novo with the 
hope that they would eventually be reduced to the already more fa- 
miliar concepts of mathematical biology. This paper represents an 
attempt at such a reduction. 

As has been discussed elsewhere (Rashevsky, 1948b, p. 564), 
reactions of an individual to any given situation may be mediated 
through a number of “parallel circuits.” Mediation through a “sub- 
cortical” circuit results in reactions which are of the nature of sim- 
ple reflexes or habits. Mediation through a “cortical” circuit results 
in reactions which are characterized by rational choices through 
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thinking. In line with the general mathematical theory of the cen- 
tral nervous system (Rashevsky, 1948b), we shall consider that all 
of these parallel circuits mutually inhibit each other. Generalizing 
the picture still further, we may consider the case in which each cir- 
cuit contains a “reinforcing” mechanism, such as described originally 
by N. Rashevsky, and subsequently developed in detail by H. .D. Lan- 
dahl (Householder and Landahl, 1945; Rashevsky, 1948b, chap. xli). 
This reinforcing mechanism results either in a facilitation or in an 
inhibition of transmission along a circuit by previous transmis- 
sions. On the basis of this picture a theory of learning has been de- 
veloped by H. D. Landahl, and successfully applied to experimental 
data. 

An individual in a society is usually confronted with the choice 
of two behavior types: One is the behavior pattern which is consid- 
ered to be “accepted” by the society; the other is any behavior pat- 
tern which goes more or less against the “secepted” standards. From 
infancy the individual is taught to behave “correctly,” that is, in ac- 
cordance with the accepted pattern. In the course of learning, cor- 
rect behavior is rewarded, wrong behavior is punished. In many in- 
dividuals this procedure results in the formation of definite habits. 
The individual acts correctly without thinking about it and without 
questioning the correctness of his actions. We may therefore tenta- 
tively assume that this habitual social behavior is mediated through 
an appropriate “subcortical” circuit. 


FIGURE 1 


Let us now consider the case in which the “incorrect” actions 
are mediated through a “cortical” circuit and are coincident with a 
rational discussion of the situation and with a logical analysis of it, 
which may lead to doubt as to the appropriateness of the socially ac- 


ie) 
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cepted “correct” reaction. This is not the only possible case, and it 
must be kept in mind that a number of socially incorrect reactions 
may be mediated through some other “subcortical” circuit, and not 
be associated with any particular reasoning. Since, however, the for- 
mer case is interesting for our problem, we shall consider it in this 
paper, without prejudice to other possible cases. 

Let us therefore assume that a “wrong” reaction not only re- 
sults in a punishment, but also puts a mechanism for rational think- 
ing into operation, as described elsewhere (Rashevsky, 1948b, chap. 
xlii), this mechanism making a rational choice of a reaction which 
may be logically preferable to the accepted one. 

Schematically the relation between the “subcortical” and the 
“cortical” circuits is shown in Figure 1, taken from our book (Ra- 
shevsky, 1948b, chap. xli). The theory of this neural net, developed 
by H. D. Landahl, and given in loc. cit., can be used here for our pur- 
poses without much change.- We presuppose the reader’s familiarity 
with the above-mentioned chapter. 

We shall, however, generalize the considerations given in loc. cit. 
by considering not only the probability of “correct”? and “wrong” 
reactions, but also the probability of “no reaction.” In other words, 
we assume that the threshold h of the reciprocal inhibition circuit 
contained in the network shown in Figure 1 is not negligible (Rashev- 
sky, 1948b, chap. xxxiv). This in its turn means that the difference 
between the stimuli S. and S, , or between the urges to act correctly or 
wrongly, must be sufficiently large in order to produce any action at 
all. Otherwise the individual remains undecided. 

By evaluating explicitly expression (7) of chapter xxxiv of loc. 
cit., we find for h > 0, and e, > Ew 


Py = fererensh, (1) 


which reduces to expression (6) of chapter xli of loc. cit. forh = 0. 

Repeating now the argument of that chapter but using expres- 
sion (1) above instead of the old one, keeping the same notations and 
putting &. — é = 0, we arrive at the following expression for the 
number w of errors as a function of the number 7 of trials 


} ] 
 k(b—8B) 


which for h = 0 reduces to equation (7) of chapter xli of loc. cit., if 
we put oc — €ow — 0 in the latter. ; 
Let «. represent the amount of excitatory factor at the connection 
8-, which corresponds to the stimulus, or urge to do as required by 
the accepted social standards without questioning. Let e denote the 


w= cag pa cet | , (2) 


og | 14 
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amount of excitatory factor at s,,. which corresponds to the intensity 
of stimulation of the “cortical” circuit, which, when set in action, 
tries out on the covert level various other behaviors which are at 
variance with the accepted one. When this circuit 7 is in operation, 
the individual does not uncritically accept the pattern taught to him, 
but subjects everything to a critical rational analysis. The threshold 
h is now a characteristic of the two parallel circuits, the “cortical” 
and the “subcortical.” The number x of trials may be considered as 
proportional to the age T of the individual. Thus, with a as a co- 
efficient, we have 


nal. (3) 
Introducing this into (2), and putting 
ka=y7, (4) 
we now obtain 
if of 
WT) iim Rpgcegy Re 1" ema er | : (5) 


Equation (5) gives the total number of errors committed by an 
individual of age T. 

As we said above, every time an error is made a set of “mental” 
trials begins until the individual reasons out a solution which is 
preferable to the accepted one (Rashevsky, 1948b, chapter xlii). An 
individual who is not capable of such a reasoning, but accepts with- 
out question the behavior taught to him, will never contribute to any 
changes in the accepted behavior patterns. A society composed of 
such individuals would never change. On the other hand, individuals 
who question the existing and accepted behavior patterns, and who 
are capable of rationally figuring out possible changes in those pat- 
terns, will eventually cause a change in at least some of them. For 
convenience we shall denote an individual who always behaves ac- 
cording to accepted patterns as conservative. An individual who ra- 
tionally behaves differently from the accepted standards we shall call 
progressive. Care must be taken not to ascribe to the above words 
their usual rather specific political connotations. 

Let the most probable number of “covert trials” required by a 
given individual to reach a rational solution of the problem be m. 
The number 1/m may be considered as a measure of the individual’s 
reasoning ability. The quantity dw/dT may be considered as a meas- 
ure of the individual’s conservatism. ; 

The “covert” or “mental” trials should not be confused with the 
“overt” trials, the number of which is defined by (3). 
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If on the occasion of each wrong reaction one “mental trial” is 
made, the most probable number of wrong reactions which is required 
to make the individual find a rational solution is m. Considering 
the rather natural case in which such an individual, once having 
found a rational solution, henceforth behaves according to that solu- 
tion, we find that the probability of an individual to begin to behave 
in this “progressive” manner at the age T is equal to 


w(T) 


m 


. (6) 


Equation (6) shows that p increases with reasoning ability. 
Equation (5) shows that w increases with T, tending asympto- 
tically to a constant, w(coc). Hence at any finite time 


w(T) <w(o). (7) 
Therefore the condition 
w(ao)<<m (8) 


will assure a very smail probability of an individual becoming spon- 
taneously progressive. 

The constant b in equation (5) represents a measure of the so- 
cial reward for a correct action, the constant f similarly represents 
a measure of the social punishment for a wrong action. 

The reward which an individual receives is in general a func- 
tion of his socio-economic status. When two individuals are behaving 
“correctly” in a society in which economic criteria play an important 
part, the one with the higher socio-economic status is likely to be 
better rewarded. Therefore it is interesting to consider the case in 
which 0 is an increasing function of the income 7, 


b=f(t). 
We may take, for example, 
b = kt or DA log (1 2). (9) 


Similarly 6 may be a function of 7. Most likely 6 decreases with 
4. For simplicity we shall limit ourselves to the case where f is con- 
stant. 

As seen from equation (5) an individual with a larger b, that 
is, in the upper income bracket, will learn the accepted behavior more 
rapidly, and also his w() will be smaller. Hence, on the average, 
individuals in the upper classes are more conservative. 

Let an individual of age T have w(T) << m, and thus prob- 
ably be naturally conservative. At that age the values of «, and éy 
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are given [Rashevsky, 1948b, chap. xli, equations (4) and (5)] by 
fo = Son UC? ff es pls (10) 
Since c + w=, we have, because of (3), 
c+w=al; col — wr. (11) 
Equations (10) and (11) give 

Ee= oc + baT —bw(T); Ep ee PUL Ye (12) 

Because of the assumption made originally that éo¢ = fw, We find 
ee — 8 = bal + (P— dy w{T): (13) 


Now consider the case in which another individual, or a group 
of individuals, tries to influence the given individual to behave con- 
trary to the accepted patterns. Formally this will have the same ef- 
fect as the addition of a positive amount ¢ to the quantity ew. 

In order for the influence to be effective, that is, in order that 
the given individual begins to behave contrary to the accepted pat- 
tern in a sufficiently large percentage of cases, we must have 


Cy On alls (14) 


Let the individual be under this influence of others for a time 7. 
The probability of a “correct” behavior, while under such influence 
is given, for éo > &, by 


Pol = 1 ekiewtp-ec), (15) 
which, if (14) holds, may be approximately written 
Pel =k elewtp-e0), (16) 
Hence the approximate probability of a wrong reaction is 
Py =1— 4 ekeuro-e0), (17) 


When ¢ is very large, or the influence is very strong, P,,' is close to 
one. 

When (14) does not hold there is an appreciable probability of 
“no reaction” (Rashevsky, 1948b, chap. xxxiv). This may be inter- 
preted as a sort of passive, inertial resistance to any social change 
without active defense of the existing situation. 

During the time 7, the individual makes, according to (3), ar 
trials, and therefore on the average 


arPy = ar [1 — § e*(euth-eo)] (18) 
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wrong reactions. If it takes m wrong reactions to reason out a so- 
cially different behavior, then, in order to make the probability of a 
permanent change to progressive behavior close to one, we must have 


ar [1 a 4 ek (Ewtp-ee) ] >™m. (19) 


From (19) we also have 
yf mM 
ees cee Les) + (S6e— 155). (20) 


Since e- — e» is given as a function of the age T by equation 
(18), equation (20) gives us the minimum influence necessary to 
“convert” an individual of age T to a progressive behavior. We see 
that ¢ increases with age. Since the expression in brackets of (19) is 
less than one, m/ar < 1. Hence expression (20) shows that ¢ in- 
creases with m and decreases with 7. 

These considerations could be generalized to the case in which 
the influence is exerted in “installments,” 7;, so that 7 = S 7;, in- 
troducing possible effects of forgetting. If the exact expression (15) 
is used instead of (16), hk will enter into inequality (20). This case 
may be of special interest. A further generalization to be considered 
is that of r different possible patterns, in line with H. D. Landahl’s 
theory of learning, to choose from any number of possibilities (Ra- 
shevsky, 1948b, chap. xli). 

Let the number of individuals at a time t whose age lies between 
T and T+ dt be N(T,t)dT. At the time t + At the number of in- 
dividuals who are of age 7 — At at t, will have an age 7. Hence, 
neglecting higher powers of A 

ON (T, t) 


Senge hee N (Ft) SUNT 1) ae (21) 


This will be the case in the absence of mortality. Actually of the 
N(T — At, t) individuals not all will reach the age 7, but only the 
fraction 

1—w(T— At) At, (22) 


or 


0 
1—w(T) at + (At)?, (23) 


where w(T) is the death rate at the age T. : 
Hence neglecting higher order infinitesimals we find instead of 


(21) 
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oN 
N(P, t+ at) =[1—w(2) 4¢ | [wer — Fat | , . (24) 


or 
oN 
N(Dyt+ At = NT, 1) ——, Ate) NATG) A (25) 
Hence 
oN 

ANN (Tt odbc WN (Di) =m, Abad TN Clete (26) 
Therefore 

oN oN 

—-=—— — w(T)N. He | 

at re tate ae 


This equation gives the variation of N(7) with time. If p(t) 
denotes the total birth rate, which as a rough approximation may be 
considered as proportional to the whole population 


N(t) = i} N(T,t)aT, 
then we have by a similar argument 
n=( ah \ (28) 
e Aged ete , 


Of the N(T)dT individuals of age T, T + dT, let the number 
6(T)N(T)dT be conservative. Each of those makes 
dw(T) 
CL 


errors per time dt, where w(T) is given by equation (5). Because 
of (6), the total loss of conservatives during the time dt is 


6(T)N(T)adT dw(T) dt 
m dt j 


dt (29) 


(30) 


This loss is due to “spontaneous” changes of individuals. 
For the probability of a change due to the influence of others we 
have from (17), putting «.— «, =A, 


adT nine 
eee Ure ph oo | , (31) 
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and therefore for the loss in the number of conservatives due to pro- 
gressive influence, we have 


a6 (T)N(T)aT ; 
pained seek eee — 1 @-k(g-A) 
m | ne | (32) 
The total number of individuals who change during the time dt 
from conservative to progressive is given by the sum of (30) and 
(32). Therefore 


_ O(T)N(T, t) dw, (T) 


—| a(T)N(T, t) | = 


ot m dT 
ab (T)N(T, t) | ey 
i {= g e-*k(g-A) A 
mM 
The quantity ¢ is given by 
g(t) = { 6(T)N(T, t)dT. (34) 


Substituting (34) into (33), we obtain a very unwieldy nonlinear 
integrodifferential equation, the discussion of which will be postponed 
until some future time. 
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Focd and metabolic waste products, insofar as they act upon the 
hereditary substrate of cells, are the most important factors governing 
tissue growth. Equations describing the growth of tissues are derived 
in consideration of this fact. A quantity is found in these equations 
which, if slightly changed, results in very great changes in the growth 
rate of the tissue, where such very great changes are interpretable as 
neoplastic growth. The relationship between our equations and similar 
equations which others have proposed is discussed. 


A tissue will be defined as a group of specialized cells differing 
in structure from other such groups in the adult body (Smith and 
Copenhaver, 1948). Our objective here is the derivation of growth 
equations for such tissues. Quantities in these equations will be 
sought which, if slightly changed, would result in very great changes 
in the number of cells in the tissue under discussion. These very 
great changes would be interpretable as neoplastic growth. The 
search for such quantities has been suggested by N. C. Wetzel (1944) 
when he mentions that abnormal growth can of itself be due to noth- 
ing more than a pathological change in one or more of the same dy- 
namic properties which govern the intricate transactions connected 
with and leading to the complete and final development of the or- 
ganism. 

Preliminarily heredity as far as it affects growth rate in the 
tissue will be regarded as a factor of constant magnitude. “H. J. Mul- 
ler (1947) states that, as a result of mitotic divisions, every cell of 
the developed individual, barring occasional accidents, will come to 
contain every gene that was present in the original fertilized egg. 
Likewise, determination will be regarded as complete in the typical 
tissue under discussion. By determination (Needham, 1942) we mean 
the fixation of the fate of a part of an organism so that—if it de- 
velops further—it will perform one kind of morphogenesis or histo- 
genesis, and one kind only, irrespective of the situation in which it 
may find itself. The values of the constants in the equations will be 
taken to represent the genetic make-up of the constituent cells as 
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well as the course of their determination. Hormonal influences will 
be adjudged, in a sense, to alter heredity, to act in the same way as 
it does, or to do both. N. C. Wetzel (1944) and M. F. Morales and F. 
L. Kreutzer (1945) suggest that the effect of hormones on growth 
may be expressed by changing the values of certain coefficients or 
constants in the applied growth equations. 

The variable by which growth will be represented is cell num- 
ber. This is used particularly since, in the type of growth known by 
the pathologist as neoplastic, the one measure of the vastly increased 
growth activity is the tremendously increased number of cells in 
some tissue. R. A. Willis (1948) defines.a tumor as “an abnormal 
mass of tissue, the growth of which exceeds and is uncoordinated 
with that of the normal tissues, and persists in the same excessive 
manner after cessation of the stimuli which evoked the change.” He 
makes it clear that what he refers to as growth implies cell multipli- 
cation by cell division. 

The mathematical treatment given below will divide tissues into 
two groups, each of which will be discussed as a separate case. Ac- 
cording to P. E. Smith and W. M. Copenhaver (1948), the duration 
of life of cells shows great variation. The nerve cells become differ- 
entiated very early and lose the power of mitosis. This is prob- 
ably true of some other cells as well, but most of the cells live a much 
shorter time. There is a constant replacement of dying and dead cells 
through mitosis of others. Such mitoses are common in the connec- 
tive tissue cells and especially in the active secretory cells. As far 
as nerve cells are concerned, they diminish in number with age. Neu- 
ronophagia occurs in the cerebral cortex of man and is found most 
abundant in the senile brain (Andrew and Cardwell, 1940). 

Case I will include tissues in which there is a constant replace- 
ment of cells throughout the life span of the organism by new cells. 
Case II will include tissues in which cells, following a short period 
after birth, no longer multiply, but undergo a gradual decrease in 
number with the advance of age. For most purposes Case I can 
readily be said to include those tissues in which cancer occurs sig- 
nificantly and Case II those tissues in which cancer does not occur 
significantly. 

Let n= number of cells in the tissue under discussion, ¢ = time, 
dB/dt = time rate of formation of new cells, and dD/dt = time rate 
of death of cells. 

Then 


Bae a a) 
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The term “medium” will be used to refer to the tissue’s immediate 
environment. This environment has been called tissue fluid (Cowdry. 
1944). Tissue fluid lies between the extravascular cells throughout 
the body and the endothelial walls of capillaries. Materials passing 
from the bloodstream must penetrate through the endothelium and 
enter the tissue fluid before they reach the cell membrane. The cells 
derive their nourishment immediately from the tissue fluid, as well 
as eliminate their wastes into it. Let F = the amount of food in the 
medium. By food we mean those chemical substances which are taken 
by the cells from their environment, whether these substances be 
solid, liquid, or gaseous. Let 7 = the amount of waste products of 
metabolism in the medium. 
In Case I we assume 


Prine Saab) 
ndt q+ (F—b) 


where a, b, ¢, and g are positive constants, and with the limitation 
that 


cT, Hieg(2) 


dB 


(3) 
dt 


y. 


; 1 
The effect of food on — 


nm aL 


to the representation of an adsorption isotherm by 


in the situation where T — 0 is analogous 


Ooo 
C0 hs eee are > 
Lis DoGy 


where o = the fraction of surface occupied by adsorbed molecules, 
¢, = the concentration of the adsorbable substance in the adjoining 
gas or solution, and b, is a constant (Hinshelwood, 1946). When Cs 
is small in equation (4), « is approximately in direct proportion to 
it, and when ¢ is large, o is almost independent of ¢. The same is 
true of the relation of (F — b) to the quantity 


a(F — b)/lq +4F —b)). 


It is to be noted that the use of b in equation (2) makes the presence 
of a certain minimum amount of F (equal to b) necessary in the me- 
dium before any cell division can occur, even when the amount of T 
present is zero. The explanation of the above-mentioned dependence 


(4) 


ROL i wal on F may not, however, be based entirely, or even signifi- 
n dt Fon 
4 cantly, on adsorption processes. N. Rashevsky (1938) shows through 
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a mathematical treatment that, while down to a certain concentra- 
tion of oxygen the rate of oxygen consumption of a cell remains con- 
stant, below the critical point the rate of consumption falls to zero 
with decreasing oxygen concentration. His discussion shows that an 
incomplete penetration of oxygen into the cell may be the cause of 
this effect. His arguments apply just as well to other diffusable sub- 
stances in the medium. The experimental work of W. Kempner 
(1936) shows that one must go down to very low oxygen tensions be- 
fore the decrease of respiration at such tensions may be perceived. 
He also showed (1937) that the growth of aerobic bacteria is in- 
hibited by lowered oxygen tension. As far as the cT term in equa- 
tion (2) is concerned, we might say that T will affect growth in low- 
er quantities than are necessary to affect other cell functions. M. 
Rubner (1924), by placing yeast cells in different concentrations of 
salt, reduced their water content and—when the concentration was 
sufficiently raised—their heat production. But only a mild degree of 
hypertony was required to stop growth. It is to be realized that the 
values for such constants as a, b, c, and q will be different in differ- 
ent tissues. Letting F get very large and letting T become insignifi- 
cantly small gives us optimum conditions for cell multiplication with 


1 dB 
are. equal to the constant a. A. Fisher (1946) says that cells pro- 


liferate at a very definite rate which is characteristic of the particu- 
lar cell type and which is in fixed ratio to the rate of proliferation 
of other types, and that differences in the duration of the interkinesis 
may be imagined to be the primary cause of organization and differ- 
entiation. Expressing the same concept in a somewhat different way, 
J. Needham (1942) mentions that a specific rate of relative growth, 
expressible by a heterauxetic constant, may be regarded as the fate 
of an organ, permitting us to speak of the determination of a growth 
rate. 


Assume 


Aas dD 
which is almost the simplest way in which de ae could be made de- 
n 


pendent upon 7. 
It is to be especially noted that through equations (5), (2), and 


1 dn 
Bae °F re is a function of F' and T only. These two factors are prob- 
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ably the two most important ones affecting growth rate. A. Fisher 
(1946) states that the primary cause of the aging of cells is not to 
be found in the cells themselves, since senility and death are due to 
external factors, i.e., accumulation of metabolic products and lack of 
nutrient substance. E. B. Wilson (1925) explains modern ideas of 
senescence to emphasize the importance of the accumulation of rela- 
tively inactive products. According to C. C. Little (1947), when 
waste material is promptly and completely removed and when a steady 
and adequate supply of food is provided, growth is continuous and 
shows no sign of slowing down. C. M. Child (1923) gives the general 
result of the aging process as a decrease in the metabolic activity of 
each unit of the organism because the proportion of the relatively 
stable constituents in the living substratum increases. J. Needham 
(1942) refers to the declining basal metabolic rate with advancing 
age as due to the changing relative proportion of respiring proto- 
plasm and inert material in the cells. 

In our present theory specific growth-promoting substances in 
the medium are not considered, except insofar as they might be ac- 
counted for as part of the constituents of fF’. A. Fisher (1946, p. 
250) says: “.... age brings about certain changes in the serum 
reflected in the rate of growth of tissue cells cultivated therein. The 
changes are characterized by an increase in the amount of growth- 
inhibiting, and not by a loss of growth-promoting factors.” The ef- 
fect of growth-promoting substances may not be important in the 
intact organism’s tissues, since these substances are produced by the 
breakdown of leucocytes, and this may not occur to any appreciable 
extent in the vicinity of most types of tissue cells. Nevertheless, it 
may well be that a more adequate mathematical treatment will ne- 
cessitate their consideration. A. Carrel and A. H. Ebeling (1923) 
mention the increased restraining effect of serum on growth in old 
age as being due, at least in part, to an increased amount of inhibit- 
ing principles. 

For our purposes it will be convenient to classify intracellular 
biochemical (or metabolic) reactions into two groups: (1) main- 
tenance reactions, which are necessary to maintain the cell proto- 
plasm in its living state; (2) growth reactions, which increase the 
amount of protoplasm in cells. Whether an organism grows or not, 
it must obtain energy for the process of living, and this energy is 
obtained from oxidative processes; some food materials are built up 
into protoplasm, while others are utilized to produce energy by their 
oxidation (Heilbrunn, 1943). J. Needham (1942) mentions that there 
are certain cases in which it has been shown that metabolism can 
occur in the absence of growth so that cell division is quite absent 
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and the cells are simply maintaining themselves as far as possible in 
the status quo. The differences between the dietary requirements of 
amino acids for growth and for the maintenance of nitrogen equilib- 
rium (Burroughs, Burroughs, and Mitchell, 1940) emphasize the 
game concept. The rate of maintenance metabolism will be assumed 
proportional to the number of living cells present in the tissue, while 
the growth metabolism can be considered proportional to the birth 
rate of new cells. In justification of this assumption, we mention 
that during the lifetime of a tissue region the cells remain within 
the same average order of magnitude as far as their size is concerned. 

We will assume that the vascular bed of the tissue region grows 
in direct proportion to the tissue region itself, whether that growth 
is of a positive or negative type. In experimental work done on 
the skin (Weiss and Frazier, 1930), it was found that the radius of 
the cell area supplied by a single surface capillary might be the same 
in an aged individual as that in a young adult, and that reduction 
in the number of capillaries cannot be held responsible for the 
changes observed with advancing age. As a consequence the ability 
of the blood vessels in a tissue to supply it with food or remove its 
wastes should be proportional to the number of cells in the tissue. 

In consideration of the foregoing discussion we set down the 
following equation: 


———=d+—-—-g, (6) 


where d, f, and g are positive constants. We let d = the amount of 
F used up for maintenance metabolism by a cell in a unit time, g = 
the amount of F supplied by the vascular bed of the tissue to a sin- 
gle cell in a unit time, and f = the amount of F used up in growth 
metabolism leading to the addition of another cell to the tissue region. 

With the exception of one additional term, the expression used 
to describe the change in 7 follows the same form as does that for 
F' , so that we have 


147 54 bdB oj aD 
oo n dt Oy EE on dt ue (7) 
where h,i,7, and m are positive constants to be defined below. The 
j dD 
oe we term accounts for waste products yielded up into the medium 


when a cell dies there. Now in equation (7) we will let h = the | 
‘amount of 7 from the maintenance metabolism of a single cell pro- 
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duced in a unit time, i = the amount of 7 produced by that quantity 
of growth metabolism leading to the addition of another cell to the 
tissue ‘region, 7 = the amount of 7 produced by the death of a single 
cell in the tissue region, and m = the amount of T removed by the: 
vascular bed of the tissue from a single cell in a unit time. 


In Case II we assume 


dB : . Hs 
dt 


at all times. From equation (1) 


—— =— — 9) 
n dt n dt ( 
Reasoning as in equations (2) and (5) 
1 dD a (Ff — 0’ 
eee 7 eee eh Yim (10) 
n at Go iO.) 


with p”, a’, q’, and b’ positive constants. In examining (10) it will 
be noted that a certain part of the death rate will be due to starva- 
tion (or lack of F'). In equation (2), which belongs to Case I, food 
deficiency lowers the birth rate, but does not affect the death rate. 
Equation (6) becomes for this case 


=== J’— gq’. (11) 


Equation (7) becomes 
—m', (12) 


with d’, g’, h’, 7’, and m' positive constants. This completes our dis- 
cussion of Case II for the present. 


In expression (2) . “ was made dependent on F’, letting what 
n 


we have called growth metabolism be a function of F. On the oe 
hand in equations (6) and (7) we have used the constants d and 

to represent maintenance metabolism independent of Ee This is es 
tified at least in Case I since in that case there is cell replacemen 
all through life indicating that there is adequate F present ig do oy 
There is considerable likelihood that, when available food is conn : 
growth metabolism may decrease and disappear though for a ‘ae 
time maintenance metabolism will continue, with a much lower leve 
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of F required to make maintenance metabolism dependent on F' than 
to make growth metabolism so dependent. Pertinent to this, A. G. 
Hogan and R. W. Pilcher (1933) showed in mammalian experiments 
that the basal metabolism of their animals was apparently the same 
regardless of the amount of protein consumed, but that the high pro- 
tein groups grew more rapidly than those on low protein rations. 

Case I is certainly the more important case for our study since 
among the tissues included within its range of description occur the 
vast majority of cancers. Summarizing Case I we list the following 
set of simultaneous differential equations. 


dn dB aD 
dt dt dt’ 
end ee pee) (ee ) 
a ae Sed Sen es 
ndt q+ (F—D) dt 
Lape 
Ree 
1 dF jf dB 
on dt nd” 
1 dT i dB j aD 
n dt n dt ndt 


It may be said parenthetically that the relation between the num- 
ber of cells in a tissue and the weight of that tissue may be given by 
the following equation 

dw dn 2 aT 

—_ = Vv —— 

ti edie ae Ce 
where w = the weight of a tissue, and where v and z are constant 


coefficients. 
Putting equations (2) and (5) into (1) we have 


1dn a(F=b) 
ndt q+(F—b) 
Using equation (2) in (6) we have 
apie nh siStat Be. 
n dt a+ (fF —d) 
Using equations (2) and (5) in (7) we have 


—cy'—*pT. (14) 


el, |=o. (15) 
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eee +i] SS er | + ior—m (16) 
n dt q+ (F—b) ; 
As is justified in our discussion about (2) we can let 
a(F —b) ; 
q+(F—v) “ a 
when F is large, but when F is small we let 
a(F — b) ECE ree (18) 
Uap li == b) q 
In this sense a large F is greater than g + b while a small F is smaller 


thang + b. 

Because of the mathematical difficulty inherent in the solution 
of equations such as (14), (15), and (16), we will use the approxi- 
mations (17) and (18). First we shall apply (17). Then from equa- 
tions (14) and (16) we have 


1 dn 
regrt ets (c+) (19) 
aT 
OA Aig Ee ee Pte Pn EP Pee TAY yh (20) 
n dt 
Let 
dx =ndt. (21) 
Then 
dn 
—=a—(c+p)T; (22) 
dx 
= h—m + ia (jp —ie)T. (23) 
x 


Equation (23) has separable variables so that we can solve for T and 
obtain | F 
es AA aes ee ag ia os m (24) 
jp —1¢ jp — te 
where c’ is an integration constant. When x = 0, let T = 0 so that 


T= depen mn | e(ip-ic)e — J] | ‘ (25) 
jp — ic 


124 TISSUE GROWTH AND CANCER 


Putting equation (25) into (22) and then solving the resulting equa- 
tion by separating variables gives 


(e+ p) (h+ta—m) 
n=| a is |# 


| jp — ie 
(26 
(c+ p)(h+tma—m) 


(jp — %¢)* 


ip—ic)x2 " 
ep vc) a Cc ; 


where c” is an integration constant. When + = 0, n = M% SO that 


(c+ p)(h + 1a—™m) (c+ p)(h +ia—m) 
Some CAO : ; a — is 
| Etna (jp — 2€)? 


(27) 
x er cecal | | =F the 3 


t= | ante (28) 
‘ n(x) 


constituting the solution of equations (19) and (20). 

Concerning the biological meaning of «, it seems best to con- 
sider it as a kind of biological time or at least a measure of such. 
The quantity x is a function of t whose derivative with respect to t 
is always positive, as is shown in its definition in equation (21), 
where n > 0 in all situations biologically meaningful. The integral 


t 
{ ndt 
0 


has been described as representing a cumulative toxic effect by V. 
Volterra and U. D’Ancona (1935), V. A. Kostitzin (1989), and C. 
N. Hinshelwood (1946). This concept is allied with A. Carrel’s 
(1931) discussion of biological time. He explains that physical time 
(= t), which is measured by a clock, differs from the time which we 
live, and that at all levels of biological organization the cause of du- 
ration seems to consist of the modification of their medium produced 
by living structures, and of the secondary changes undergone by these 
structures under the influence of the modified medium. The same sort 
of idea is emphasized in.a different manner by F. G. Donnan (1936), 
when he mentions that for the description of the change of living 
systems in time a pure differential analysis will not suffice. In such 
cases one must use instead integro-differential equations giving a 
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Ee . . mathematical description of a living organism as a special 
type of world-tube in space-time.” 

- The quantity jp — ic is the difference between two quantities, 
which are surmised to be of the same order of magnitude. If this 
supposition is correct, then small changes in jp or tc, or-in both, 
which will permit each to remain of the same order of magnitude as 
it was previously, will cause relatively much larger changes in the 
magnitude of jp — ic relative to its original magnitude. Attention 
is now directed to determining how changes in jp — ic could account 
for a change of growth from the normal to the neoplastic type. It 
should be remarked that we certainly do not believe that changes in 
j,p,%t, or c are the only ones which could account for a transforma- 
tion to cancerous growth. Although we do not engage in a detailed 
study of them at this time, positive changes in @ would seem quite 
capable of producing neoplastic proliferation. 

-Taking the derivative of-equation (27) we obtain 
ee iat ss (c+ p)(h+ia—m) | Pee wanes 


: : (29) 
da 4) — te 


Setting equation (29) equal to zero and solving for %mae, it is found 
that 


Xmas 


BUDE) +1]. 
(c'+ p) (h+ta—m) 

CF h+ia—m) ‘ 

If 7p — ic is less than — Cee Lmaz Goes not exist. 
In normal growth there is a positive %mox, Which gives a maximum 
nm. Then normally 7p — ic must be greater than 


be (c+) (h+1a—m) 
a 


The condition for the existence of a maximum of is not only 
that dn/dx = 0 be an equation having a positive root x = Lmaz » but 
also that d?n/dx? be negative for * = Xn. Taking the derivative 
of equation (29), we obtain 


d?n/dx? =— (c + p) (h'+ ia — m) e492 , 


In order for d?n/dx? in the above equation to be negative for any 
value of x, (c + p) (h'+ ia — m) must be positive. Since both c and 
» must be positive from their biological meanings, then, from the 
above-mentioned considerations, (h + ia — m) must be a positive 
quantity in order that normal growth may occur. 


ser 
jp — 
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We have carried out a graphical integration of the expression 
whose integral is indicated in equation (28). In order to do this 1/n 


. . 


was plotted against x in Figure 1 with jp — tc given different values 


Yn, 


FIGURE 1 


in each of the figures. (In all the figures the symbol y is to be taken 
as having a meaning identical with the expression jp — ic.) In Fig- 
ure 1, when 1/n was plotted against « when ip — ic=0, it is noted 
that an indeterminate expression of the form 0/0 had to be evaluated, 
as is obvious from the form of the functional relationship between 
jp — ic and n given in equation (27 ). Returning to (23), letting 
jp — ic = 0, integrating the resulting expression, and assuming as 
before that when « = 0, T = 0, we have 
T=(h—m+ia)z. 
Substituting this result into equation (22) gives 
dn/dx=a—(c+p)(h+ta—m)z, 

which when integrated with respect to ~ becomes 


n=ax— (c+ p)(h+ 1a— mM) x?/2 + N. 


y=0.02 


\ 2 
FIGURE 2 
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This is the equation we have used in plotting th i i 

The planimeter was utilized to ee nese ees ee 
the graphs enabling us to plot ¢ against x in Figure 2. By using tiene 
graphs (Figures 1 and 2), it is possible to find values of n and t 
corresponding to a given x so that a direct relation between n and t 
is obtained, with the variable x eliminated. This relation is graphed 


FIGURE 3 


on Figure 3. In preparing the graphs arbitrary values have been as- 
sumed for a and (c + p)(h + ta — m). (In preparing the graphs 
given here we have used n= 0.5,a = 38, and (c + p)(h + ta —™m) 
= 1.) To obtain actual values for a, (c + p)(h + ia — m), and 
jg — ic on the basis of available data requires a tedious numerical 
computation which we have not done at the present time. It should 
be remembered that the units used in the figures are purely arbitrary. 
From Figure 3 it can be seen that, when 


jp —ic >— (c+p)(h+ia—m)/4, 


the resultant curve giving n as a function of t is quite representa- 
tive of the pattern which normal growth takes. This curve fits well 
the description of normal growth given by D. W. Thompson (1942). 
He states that rate of growth, which is represented by the slope of 
our curves in Figure 3, depends on the age of the organism. It has a 
maximum early in life, after which epoch of maximum it slowly de- 
clines. This means that there is a point of inflection, or abrupt change 
from an increasing to a decreasing velocity of growth. Negative 
growth occurs with old age. The marked manner in which growth is 
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affected, in the way of acceleration or retardation, at certain physio- 
logical epochs of life, such as birth and puberty, is explicable on the 
basis of a change in the values of the coefficients due to hormonal in- 
fluences at the time of puberty and to completion of determination at 
the period shortly following birth. 

It is seen in Figure 3 that a decrease in jp — ic produces an in- 
crease in growth rate. If jp — ic is decreased down to 


—(c+p) (h+ta—m)/a, 


n no longer has a maximum, but approaches an asymptote represent- 
ed by the line n = a?/[(¢ + p)(h + ia —m)] +7. Though we have 
not plotted a curve for jp — i¢ ae Cee) (h + ia— m)/@,a study 
of equation (29) shows that when much time has elapsed, the e/?*°? 
term is insignificant and ” follows: the so-called exponential type of 
growth, in which dn/dt is equal to some positive constant times ”. 
This curve never has a point of inflection at all. 

A small decrease of jp — ic from the normal would cause the 
new growth curve to reach a maximum larger and later than, but 
close to, the normal maximum. When the decrease in jp — ic is great, 
the very large maximum » may never be reached during the life span 
of the individual, but occurs at the end of a longer time; or, if the 
decrease in jp — ic is sufficiently great, there will never be a maxi- 
mum 7. The former type of growth may be interpreted as that of 
a benign tumor and the latter as that of a malignant tumor. 

We have already made some mention of the factors determining 
the magnitudes of our constants. Of these, hormonal and genetic in- 
fluences are variable in the postnatal lifetime of an organism. As 
far as endocrinological activity, this is obvious, but genetic factors in 
somatic cells undoubtedly change as well. Important in this is the influ- 
ence of stray radiation and the molecular readjustments occurring in- 
cident to the thermal agitation present at the temperatures at which 
life exists (Schrodinger, 1945; Muller, 1947 ). The concept ennun- 
ciated here assumes the existence of somatic mutations, which give 
rise to changes in the individual body which, of course, cannot be 
transmitted by the germ cells, but may be perpetuated by vegetative 
reproduction (Woodruff, 1941). Also causing alterations in the mag- 
nitudes of the constants very early in life are the changes leading to 
the completion of determination. That determination is not complete 
until some time early in postnatal existence is suggested by the in- 
ability of the newly born organism to form antibodies and the corre- 
spondence of this inability to the time when antigenic specificity 
among tissues has not yet developed (Needham, 1942). This situa- 
tion is significant with regard to the causation of embryonic tumors, 
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in which the disturbed chemistry of the organizers with ‘its conse- 
quent effect on the course of determination is most important (Wil- 
lis, 1948). = 

When a change in jp — ic occurred, it might be argued ‘that it 
would take place in only a very few cells in a tissue, while the: pres- 
ent theory has been described as giving the behavior of a tissue as a 
whole, with n being the total number of cells in the tissue. In order 
to apply the theory to this situation, we might regard jp —icas an 
average over the entire tissue, with any change occurring: gradually 
and reaching its final value when the entire tissue is replaced by the 
cells possessing the changed value of jp — ic. An alternative would 
be to use a different ” in this situation, where the » would:not be the 
total number of cells in the particular tissue, but would be the num- 
ber of those cells with the altered value for 7p — ic. Here again, as 
the tissues were replaced by neoplastic (having the lowered value of 
jp — ic) cells, this n would become identical with the n described in 
the general definition of n mentioned earlier. r 

Without further study it might seem obscure just how the pre- 
viously mentioned changes in such constants as 7, »,7, and ¢ could 
result in sudden increases in growth rates. Remembering equations 
(2), (5), (7), (17), (21), and (25) we obtain 


dT /dz = (h + ia— m) er-ioz ; (30) 
dD/dx =p[(h + ta —m) / (jp — ic) ] [e749 —1];,. (31) 
dB/dx =a—c[(h + ia—m)/ (jp — ic) ] Leo? 19? — 1}. zc (32) 


No matter how big or how small jp — ic may be, dT/dt is always 
positive. A positive jp — ic will cause d?T/dx? to be positive, while 
a negative jp — ic will cause d?T/dx* to be negative, with aT /dx 
approaching zero as t increases. No matter what the magnitude of 
ip — tc may be, dD/dt is always positive. When jp — ic > 0, 
d?D/dxz? > 0; when jp — ic < 0, @D/dx? < 0 with adD/d« approach- 
ing a constant value — p(h + ta — m)/(7p — ic) as t increases. 
When jp — ic > 0, dB/dt can become negative when ¢ is great enough. 
This is prohibited by the limitation imposed in equation | (3). In or- 
der to take care of this difficulty, we must put a limit on the size 
which a positive jp — ic may possess. Considering the fact that x 
reaches a maximum when v has come down to a ZeLO value at some 
late time and remembering equation (27), it is evident that’ this 
maximum permissible value for jp — ic is then the proper positive 
solution of the following simultaneous set for (7p — 40) on 


el (h + ia—m)/(ip— ie) nJ [erste —1] =a, (88) 
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[a ts (¢ +p) (h +4ia— m) / (jp iP. 16) mJ Lm 

—[(c +p) (h + ia —m) / (Jp — 26) m7] [etir-te) mem — 1] +2, —0 een. 
Equation (32) further shows that the inhibition to cell birth (or cell 
division) per cell by T approaches a constant limiting value c(h + ta 
— m)/(jp — ic) when jp — ic < O ast gets greater. An examina- 
tion of equation (32) shows that there is a minimum value down to 
which jp — ic can go without making dB/dt < 0, depending on the 
maximum size which x can eventually reach. This size of x is a SO- 
lution of equation (34) for x, . This minimum value for jp — te will 
be obtained as a negative root of (33) when it is solved after the 
proper solution of (34) for 7» has been substituted into it. 


Using equations (29) and (1) we have 


dn/dx = dB/dx — dD/da 
=a—[(e +p) (h + ia—m)/ (jp —ée)] erro" — 1) 


=a+ (¢ +p) (h + ia—m)/ (ip —%) ae 


= [Cc + p) (h al ia —m) / (79 —ic)] [ertoe] : 


From the above manner of expressing dn/dx , it can be seen why, when 
jp — ic becomes less than a certain minimum value — (c+ p)(h+u 
—m)/a, dn/dt never becomes negative. When x is large, — [(e + p) 
(hit ta —.m) / (9p — ic)] e%»-#2 becomes insignificant, provided 
jp — ic is negative. If jp — ic > — (c + p)(h + ta — m)/a but is 
still negative, then a + (c + p)(h + ia — m)/(jp — ic) < 0 and so 
causes dn/dt to become negative when ¢ is great. However, when 
ip —ie <— (c+ p)(h + ia — m)/a,a + (¢ + p)(h + ta —m)/ 
(jp — ic) is positive and so when t is great and — [(¢ + p)(h + ta 
— m)/(jp — tc)] e%?-4o* is very small, dn/dt is still positive. It is, 
of course, quite evident that as long as jp —ic <0,—[(e +p) (h+ 
ia — m)/(jp — ic)] e”*>* is a positive quantity. 

Taking the derivative of equation (19) with respect to x gives 


d2n/da?=— (e+ p)dT/dx. (36) 
) 
Using equation (23) in (36) we have 
d'n/dic? = — (c'+ p)[h+ia—m + (jp —ic)T]. (37) 


Solving equation (19) for T and putting this expression into (37), 
we get 
‘dn/dxu?=— (c+ p)[h—m + ta 


+ (jp —ic) (a—dn/dx) /(e + p)]. wr 
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Integrating (38) with respect to x and remembering that when T=0, 
dn/dx =a, and using (21), we obtain the following: — 


dn/dt = [a—n, (jp — ic)] n— ‘Gratiaye : 


e (39) 
—((e +p) (h—m + ia) + a(ip—ie)}n [nat 
When : 
jp —ic=— (0+ p) (h—m-+ ia) /a,. i, (40) 
dn/dt = [a— n,. (jp — ic) ] n— (ie — jp)n?.\> ° (41) 


Equation (41) represents the logistic curve. A large amount of sta- 
tistical and experimental study has shown that populations of the 
most diverse kinds are found statistically to follow in their growth a 
particular type of curve, the logistic (Pearl, 1927). 

A rationalization of the logistic curve was made by A. J. Lotka 
(1925) in the following way. He assumed that 


aX /dt— F(X), (42) 
where X was the mass of the system. The equation F(X) — 0 has a 
root at X = 0, for at least one individual is required to start the 


growth of a population. Expanding by MacLaurin’s formula, there 
will be a series for F lacking the constant term. Thus 


dX/dt = F(X) =aX + 0X? + oX3----, (43) 


Furthermore, equation (43) should have at least one other, root, since 
there must be some upper limit to the growth of the population. The 
simplest case satisfying this condition is that in which the right-hand 
member of (43) terminates at the second degree term, i.e., 


dX /dt =aX + bx?. . + (44) 
V. A. Kostitzin (1939) derives the logistic curve from somewhat 
different considerations. He assumes thatt | 
Number of births in (t,t + h*) = Np*h’, 
Number of deaths in (t,t + h*) = Mp‘h’, 
where he calls N and M uncorrected coefficients of birth rate and 
mortality, and where p* is the population. From this we get 
p(t +h) —p*(t) =Np*h* — Mp*h’. 
Dividing by h* and proceeding to the limit as h* approaches zero, we 
get 


+Asterisks are used after all of Kostitzin’s symbols which might be con- 
fused with our symbols. 
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p* (t) =Np*(t) — Mp*(¢). (45) 


The coefficients N and M are not constant but are dependent on p*. 
Kostitzin puts this dependence in the following form 


N=n*— vp", (46) 
M==m* c2np%, (47) 


where n‘, v, m*, and u are constants. Putting equations (46) and 
(47) into (45) gives us 


p* (t) = ep*(t) — hp" (¢), (48) 
where © 
. e=n*—m"*, (49) 
and 
h=yv+u. (50) 


In order.to show the equivalence of Kostitzin’s work and our 
own, it is. necessary to substitute the following for equations (46) 
and (47)... 

NES oy (PF Pays (51) 
M=m* + u(p" — po) ; (52) 


where'p,; is the value of p* when t = 0. Using equations (32) and 
(51) ‘it becomes evident that by substituting t = 0, we obtain 


Gwe 


By using equations (31) and (51) and letting t = 0, it can be seen 
that 


m0". 

From equations (41), (48), and (50) it follows that 
*=1¢e—jp=vt+ up. 

If p mm 0 , adD/dax = 0 and u=0. Then we have shown that 

ele. y=1C. 


Consequently, when p has a value greater than zero, it must be true 
that 


u=—jp. 
Equation (39) is of exactly the same form as the equation pro- 
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posed by V. Volterra and U. D’Ancona( 1935) and solved exactly by 
V. A. Kostitzin (1939) to give a solution having the same mathemati- 
cal form as the one we have obtained in (27). The 


t 
if ndt 


in equation (39) is described by Kostitzin as representing a total 
toxic action on the birth and death rates because of the accumula- 
tion of metabolic products. We can place Kostitzin’s equations (46) 
and (47) into our own notation and also show in the more complete 
case that 


t 
N=dB/dz =a + ien,—ien —c(h—m) f ndt , 
0 
t 
M = dD/dx = jpn, — jpn + p(h—m + ia + aj ) f Nate 
0 


The means by which the coefficients for the 


t 
{ ndt 
0 


term were selected become obvious if we assume a zero death rate 
due to » = 0 in equation (39) and then obtain the effect of the 


t 
f ndt 
term on the birth rate. 
V. Volterra (1938) proposes the equation; 


adN/dt=— (c'—b*°N + @*N?)N. bis (53) 
If we assume as an approximation that 
dt=rdn, 

where r is a constant, and use this approximation to carry out the 
indicated integration in (39), we have 


dn/dt = [a—n,(jp—ic)] n— (te — Jp) n* 
—{(c+p) (h—m + ia) + a(jp —ic)} (mM —n,°) /2, 


which is of the same form as Volterra’s equation. 
In a previous»paper (Kesselman, 1948), we proposed a growth 
equation. The G; used there is equivalent to a in the present treat- 


+Asterisks are used after all of Volterra’s symbols which might be.confused 
with our symbols. 
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ment, while the J; is the equivalent of nT(c + p). The discussion of 
the nature of the action of carcinogenic agencies given there, which 
stresses the importance of the relative magnitudes of NiG; and Ji, is 
just as applicable to the present theory. 

If we assume that neoplastic growth signifies an increase of dn/dt 
above the normal value, Figure 3 shows that two tissues of the same 
histologic type with different values of y may grow at the same rate 
during early life, but later one of them may change to neoplastic 
growth, the other maintaining the normal growth pattern. Changes 
in y occurring later in life must be examined in a different way, since 
growth curves of tissues with different values of y are then quite in- 
dividual in character. 


n 
A wm Y=-0.33= Vy 
7 Wa, =-0.33=Y. 
2) my 
RO A 
6 te x 
. Le Wa, 
ag Y= -0.33=Yy, 
e A 
x 
Ys w y= -010=Yy 
3 
m= 0.00 = Yay 
2 I Ir y= 0.02 =Yoy 
| ¥=0:50 = ¥. 
z y=0J0=¥, 
t 


{ 
FIGURE 4 


Suppose that at time t, a change-in y (see Figure 4) from yw 
to Yvr (given as Yur, in Figure 4) occurs in a tissue composed of 14 


cells at £,. Curve VI will not describe the growth of the altered tis- 
sue after t,, since at t4 tissue VI (described by curve VI) will con- 
tain a number of cells much greater than n,. Let us cut away at t, a 
portion of tissue VI so that the cells remaining in it are m4 in num- 
ber. At this time it becomes necessary to realize that in the present 
study there has been a tacit assumption that the volume occupied by 
T and F has remained constant in time. This is obvious on examina- 


tion of equations (2) and (5). In these equations neither 2 oe or 


1D te 
rs can depend on total quantities of T or F,, but can only be de- 
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pendent upon concentrations of these materials. Suppose that n = ny; 
in equation (14) giving 

os dmy,; a(F—b) 
Nyt at qd ae (F— b) 


eT — pT. 


If 1/s of ny; is removed from tissue VI, then 
1 (7 5) Sool =p) 


—cT — pT | 
Ny1/s , dt q+ (F—b) . ye a) 


since this partitioning will not affect the right-hand side of equation 
(54) at all. It becomes obvious that tissue VI, has its growth after 
t, correctly described by 


poss é ested) Ks 
ine { i= (F —d) (c+ p)T | amd, 


0 


where s is the ratio of n, to the number of cells in tissue VI at t,, 
and where the integral expression has already been plotted as curve 
VI. In Figure 4, Nyt , has been plotted against ¢. 


From the above discussion and an examination of Figure 4, it 
appears that a decrease of y in early life will not alter dn/dt much, 
but in later life will have a far more marked effect in doing this. 
This observation agrees with the observed steep and continuous rise 
of the incidence of fatal carcinoma in old age (Willis, 1948). The 
graphs further suggest, in agreement with certain experimental and 
statistical observations, that cancer is a disease of the elderly, be- 
cause of the unusually long Jatent periods elapsing between the ap- 
plication of carcinogenic stimuli and the development of tumors. 

Tumors in early life tend to progress more rapidly and dissemi- 
nate more frequently than do the same tumors in later life (Moore, 
1944). This means that in a young person a cancer grows much more 
rapidly than in an older person. Suppose that tissue IV at ¢ = tz un- 
dergoes a change in its value of y from Yr to yy:. The altered growth 
curve of this tissue has been plotted as curve VI; in Figure 4. Now it 
is possible to determine what the growth rate of a neoplasm will be 
if the decrease in y occurs early in life, when t = tz , and when t= ft, . 
Looking at Figure 4 makes it apparent that for the same carcino- 
genic stimulus (the same decrease in y) there is a greater rate of 
growth of the neoplasm early in life than there is later in life. 

In the following we shall use the approximation (18) in carry- 
ing out the solution of equations (14), (15), and (16), assuming 
that F is small. Then it is found that a 
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dn/de=a(F — b)/4a— (e+ p)T; (55) 
—dF/dx =d—g + fa(F —b)/q—eiT; (56) 
dT /dx=h—m + ia(F —b)/q + (jp —ic)T. (57) 


It should be noted in this set of equations that when the solution for 
F becomes greater than gq + 0, we should then apply the set of equa- 
tions (19) and (20). 


Let 
a/q=a,—ab/q=y,— (e+ P) =p6,g—d+afb/a=¢, 
of =o, h—m—iab/q=y, ia/g=1,ip—te=s, (58) 
__4f 
p ares, = 
Then 
dn/dx=oF + pT + y; (59) 
dF /dx = pF + oT + 4; (60) 
aT /de=iF + 6T iy. (61) 
Solving equations (60) and (61) simultaneously it is found that 
PT /dx? + a,dT/dx + T= a3, (62) 
where 
— (p+ 6) =a, , dp — da = Ay , Ap — pp = as - (63) 
The general solution of equation (62) is 
T= Oy Oped OO aCe es (64) 
where 
m, = (— 4, + Va, — 42) /2, 
and 


Ms — (—a,— \/ a,?— 4a) /2 ’ 
and where c, and ¢» are integration constants. Solving equation (61) 
for F and substituting the value of T given in (64) into it we have 


+ em (CyMz — C26) /d. cee) 


Using equations (64) and (65) in (59) we obtain an expression of 
the form i 
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dn/dx = k, + k,e™* + k,em* , (66) 


where k,, k., and k; are constant coefficients. Integrating (66) we 
find that 


N= ky + kya + kpe™*/m, + k,e™*/m, (67) 


with k, an integration constant. For the present we defer the study 
of the properties of equation (71). 

V. A. Kostitzin (1939) objects to the use of a single equation to 
represent growth, holding that the process must be divided into a 
number of stages, each characterized by its own values of the coeffici- 
ents in the growth equation used. We have given our coefficients 
variable magnitudes in a number of cases. The solution of such dif- 
ficulties is to be accomplished by the use of differential equations in 
which the coefficients are considered variable. Another possible ob- 
jection to the present theory is that it does not take into account in- 
crease in average cell size leading toward the direction of maturity. 

The author is deeply indebted to Dr. Anatol Rapoport for stimu- 
lating discussions and guidance in helping him to clarify his ideas on 
many points; to Dr. I. Opatowski for a critical reading of the manu- 
script and his helpful suggestions; and also to Professor N. Rashevsky 
for his continued inspiration. 
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Analysis is made to show that the mitotic index is simply proportion- 
al to the ratio of the duration of mitosis (T) to the intermitotic time only 
under special conditions. In the case of exponential growth of cell popu- 
lation the simple proportionality holds if the product of T and the 
growth constant is small. For power law (t”) growth of cell popula- 
tion the simple proportionality holds only when a steady state of growth 
has existed for at least ten intermitotie periods. The simple propor- 
tionality does not apply in conditions of transient growth. 


The mitotic index, m/N , for tissues in which mitosis is occurring 
is defined as the ratio of the number of cells, m, seen in mitosis to 
the total number of cells, N , examined. The mitotic index has been 
used to determine the duration in time of the process of division 
(Hoffman, 1947; Gray, 1931; Laughlin, 1919; Wright, 1925) on the 
assumption that it is directly proportional to the ratio of the dura- 
tion of mitosis to the total time between mitoses. In the following 
analysis we shall show that the direct proportionality holds true only 
under special conditions. 

Designating the mitosis duration as T and the total intermitotic . 
time as L, then m/N = T/L when the value of N is constant with 
time, i.e., as many cells are being removed as are being produced. 
Contrasted with this case of the stationary population there are the 
cases in which the cell numbers are known to increase exponentially 
with time as in transplantable mouse tumors (Hoffman et al, 1943; 
Reinhard et al, 1945) and in regenerating liver (Brues and Marble, 
1937). The purpose of the following analysis is to show how the 
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mitotic index depends on the parameters of the life cycle (Figure 1) 
of the cell and on the rate of growth of the cell population. Case I 
is concerned with the tissue in which exponential growth of cell num- 
bers occurs, while Case II deals with the general case in which the 
cell number increases as the nth power of the time. 


MITOSIS TIME 


INTERMITOTIC PERIOD >! 


t=O t=P t=L 
EARLY OR, t=T+P 
PROPHASE TELOPHASE 


COMPLETED 


Ficure 1. Parameters in the life cycle of a single cell which undergoes di- 
vision. The time between the complete formation of daughter cells is L. When 
T is small, L is approximately equal to P, the intermitotic time. 


Case I. The number of cells in the tissue mass increases expo- 
nentially with time. 


Suppose that no cells are removed. Examples of this case are 
given in Table I. The number of cells present is then N,e* where t 
is the time in days, N, the initial number of cells present, and 4 is the 
characteristic exponential constant. 


TABLE I 
Tissues in which Cell Numbers Increase Exponentially with Time 
TYPE OF TISSUE nN REFERENCE 
Transplantable mouse adenocarcinoma 0.37 days-1 geez 
Ascite Tumor cells in mice 0.6)” 3 
Transplantable mouse tumor Oe? 4 
Regenerating liver tissue (rats) eo 2) Re 5 


*Maximum, initial value of A. 


J. G. Hoffman et al, 1943. 

M. C. Reinhard et al, 1945. 

H. Lettre, 1948. 

S. Konsuloff, 1933. 

A.M. Brues and B. B. Marble, 1937. 


eS BS GOS p= 


Table I shows the order of magnitude of the values of 4 as meas- 
ured by the volume rate of growth. It is assumed that the cell vol- 
umes of the tissues concerned are restricted to fixed limits. There- 
fore, the value of 4 applies to the growth of cell numbers as well as 
of the tissue volume. . 

The analysis of the mitotic index is carried out as follows: Sup- 
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pose that the time schedule in the life of a tissue cell is given by Fig- 
ure 1. The newly formed cell grows for a time, P, then undergoes 
mitosis for a time, 7; the total time between the appearance of the 
daughter cells being L. At a time, t, let there be N cells present. 
The number, m , of cells in mitosis has an age between P and P+ T. 
All of these cells must have been born at a time between t—L and 
t—P or t-L'+ T. The rate of births at any time is 2(dN /dt) for 
binary fission. Therefore, the value of m is given by: 


= 2N (ee?) — e(t1)\ . 
or (1) 
1 ( 2er(t-L) NV, (er? oa, 1) é 
The mitotic index is m/N and thus becomes: 
m/N = 2e"" (7 — 1). 


But, e** = 1/2, since L is the number doubling time for the cells. 
Therefore 


m/N = (e&? —1). (2) 
Since AT is usually small, this can be approximated by 
m/N A271 —*since. 6 = 1+ AT , (3) 


and if L is the doubling time, the growth constant, 1, is related to L 
by: 
jin 2 Ls 
and the mitotic index becomes 
m/N = 0.693 T/L. (4) 


The essential feature in the expressions (3) and (4) is that the index 
m/N is proportional to the ratio of the mitosis time, T , and the dou- 
bling time, ZL. And in this case, the value of 4 is a measure of the 
volume rate of growth as well as of the rate of production of cell 
numbers. Equation (2) shows that if 27 is large the mitotic index 
is no longer simply related to the ratio T/L. However, equation (2) 
permits an estimation of T once m/N and 4 are known. 


Case Il: The number of cells in the population increases as ee 


The analysis of the case is carried out similar to that given in 
Case I. The number of cells is given by 
NN lat +1) 9eis (5) 
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the doubling time, t. , is given by NV /No = 2, or 

foo) 1/0 (6) 
The mitotic index for the general case can be shown to be: 


m (t—b +T7)»— (t—L)” 


(7) 
N PCa 


By expanding both terms of the numerator and the second term of 
the denominator it is easily seen that for very large values of t we 
have 


m/N=aT. (8) 


Discussion. 


An outstanding feature in the foregoing discussion of the mito- 
tic index is the ratio T/L. Few data exist concerning the duration 
of mitosis, 7, which is the datum required if the mitotic index is to 
be used for determining L. Turning to quantitative data on expo- 
nential growth of a tissue mass we shall estimate the duration of 
mitosis, 7, for the Ascite tumor cells in mice from the mitosis data 
as reported by H. Lettre (1948). He gives values of the cell number 
doubling time, L, as 25.7, 27.7 and 28.8 hours. From his Table I, 
we find three values of the mitotic index: 1.5, 1.0, and 1.1%. To com- 
pute 7 we shall use the mitotic index or (mi) = 1% and L = 25.7 
hours; T = (mi) X L/0.693 = .01 X 25.7/0.693 = 0.4 hour or 24 
min. This is the lowest value of T to be derived from Lettre’s data. 
Lettre estimates that 7 = 15 min. from the rate of build-up of mito- 
ses in the tumor due to colchicine poisoning. 

The transplantable mammary adenocarcinoma in mice described 
by J. G. Hoffman et al (1943) and M. C. Reinhard et al (1945) lends 
itself to quantitative estimates for the mitosis time on the basis of 
the mitotic index. The values found (unpublished data) for the in- 
dex were: (mi) = 0.3 and 0.5%. The growth constant (Hoffman, et 
al, 1943) was 4 = 0.37 per day. Hence, T = 0.004/0.37 = 0.0108 days, 
or 15.5 minutes, where we have taken 0.4% as the average value of 
the mitotic index. 

A. M. Brues and B. B. Marble (1937) in their study of the ex- 
ponential growth of regenerating liver estimate an average time du- 
ration of mitosis of 49 minutes. On the other hand, in human scalp 
epidermis, where it is reasonably assumed that growth is not expo- 
nential, J. M. Thuringer (1928) estimated the mitosis duration as be- 
ing from 15 to 30 minutes. : 
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The analyses given in Cases I and II indicate the nature of the 
variables which may enter into an interpretation of the mitotic index. 
For example, in early squamous cell carcinoma of the cervix or of 
epidermis, where the mitotic index averages 4% and is sometimes as 
high as 10%, it is first necessary to determine the possibility of T 
being very long (greater than 1 hour). If T is very long then the 
mitotic index is correspondingly increased and does not necessarily 
indicate a high rate of growth. Next it is necessary to know if any 
of the cells in the preparation were near a region in which they could 
be removed by differentiation into another cell form or by direct lysis. 
If neither of these possibilities exists it may be assumed that the 
population was growing exponentially. If there is evidence that cells 
might be removed by these processes it is necessary to estimate by 
which power law of the time, t, the cell population is increasing (or 
indeed decreasing, depending on the tissue region seen). A similar 
problem in estimating cell removal rates occurs in the use of the mi- 
totic index to measure growth in tissues under hormonal stimulation 
such as in the proliferation induced by estrogens. 


Summary: 


1. It-is shown that for exponential growth of a cell population the 
mitotic index is simply proportional to the ratio of the duration 
time of mitosis to the intermitotic time only when iT is small. 
For power law (¢") growth of a cell population the mitotic index 
is simply proportional to the ratio of the duration time of mitosis 
to the intermitotic time only when the population has grown for 
times long compared with the intermitotic time. 


2. For exponential growth of tumors the experimental data are used 
to compute the time duration of mitosis. 


8. Special cases in which the mitotic index cannot be applied due to 
transient conditions of growth stimulation are discussed. 
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BOOK REVIEW 


' NORBERT WIENER. Cybernetics or Control and Communication ic RhepArinal and 
The Machine. 1948. 194 pp. New York: John Wiley & Sons, Inc. Paris: Her- 
mann et Cie. a hee 


Cybernetics is the name coined by Professor Wiener and ‘his ‘friend, Dr. 
Arturo Rosenblueth, for “the set of problems centering about ‘communication, 
control and statistical mechanics” (p. 19) and the present “preliminary book” 
(p. 153) gives the fundamental concepts on which this branch: of mathematics is 
based as well as a brief survey of the fields to which it might: be applied. 

The introduction of about 30 pages relates in detail how cybernetics grew 
up. It tells how Wiener, Bigelow and Rosenblueth came to publish: ‘their classi- 
cal paper, of the Princeton meeting of physicists and physiologists; how, then, 
the Macy Foundation made possible regular meetings in which the Princeton 
group was enlarged by psychologists, sociologists and anthropologists as well 
as by neuroanatomists and neurophysiologists. It tells how the author was re- 
quested to write the present book, and it hints finally at the sociological conse- 
quences when more and more controls will be built by applied' cybernetics which 
are “bound to devalue the human brain at least in its simpler and more routine 
decisions” (p. 87). ae 

The word “Cybernetics” is derived from the Greek « v Be p.0 9:7 7 o the 
latin gubernator. The present reviewer regrets the rendering of the Greek “K”, 
the Latin “G” by a “C”. He would have preferred andther: K.Of course, there 
are cathodes, but a patient suffers from a kyphosis! © ...°! (37: 

The first four chapters of the book explain on 96 pages the theory. It be- 
gins easily enough, with a charming quotation from WilHelm:Hey (d...1554) : 


“Weisst Du wieviel Sternlein stehen ety 
An dem blauen Himmelszelt? | ae oy, - 
Weisst Du wieviel Wolken gehen " 
Weithin tiber alle Welt? 

Gott der Herr hat sie gezahlet 
Dass ihm auch nicht eines fehlet =e 
An der ganzen grossen Zahl.” ee eee, 


This leads easily to the recognition that the meteorologist‘deals ‘with a vast 
number of approximately equal particles, closely coupled ‘together,' while the 
astronomer studies a relatively small number of particles, greatly diverse in 
size but coupled with one another only in a very loose’ way. “The ehapter goes 

“on to discuss the direction of time, evolution, automatons and! robots through 
the ages, all in a simple language which even a biologist can understand. 

The steep ascent te the pinnacle of mathematics is undertaken’ in’ the next 
chapters. The second one, starting out from Gibbs’ statistieak nvechanics and 
Lebesgue’s theory of measure shows “how the whole theory of Creeus Uae RESET 
ing transformation can be reduced to the theory of ergodic ‘transformation. 

_ The third chapter defines information on’ the premise’ that “Gniformation is the 


145 ak 


146 BOOK REVIEW 


recording of a choice between two equally probable simple alternatives.” The 
amount of information can be defined as “the negative logarithm of a quantity 
which we may consider a probability” and “is essentially a negative entropy” 
(p. 78). The theory of time series, information and prediction is then developed 
at some length. The fourth chapter concerns itself with the theory of feed-back 
mechanisms, their role in control of performance, and the consequences of their 
breaking down into oscillations. 

The rést’ of the book shows how cybernetics is applicable to biology and 
sociology: Computing machines and the Nervous System, Gestalt and Universals, 
psychopathology and, finally, language and society are examined from a cyber- 
netical point of view. Throughout it is stressed that the important problems are 
not those of energy. or of matter but those of order, i.e., of information. “The 
mechanical brain does not secrete thought as the liver does bile as the earlier 
materialists claimed, nor does it put it out in the form of energy, as the muscle 
puts out its activity. Information is information, not matter or energy. No 
materialism which does not admit this can survive at the present day.” 

This truth may. by now be self-evident to many, if not to all, who actively 
work on problems concerning the brain, and is implicit in the whole work of 
Rashevsky. and his school as the numerous diagrams of neuronal circuits which 
they have constructed show. But the tendency to chase “nervous energy,” “Sin- 
nesenergie” through. the labyrinths of cerebral pathways is perhaps still lin- 
gering on among: epistemologists and their friends. Wiener’s book can be ex- 
pected to clarify.definitely our thoughts. 

Space does not permit to discuss the last chapters in detail. There is much 
food for thought in them, and there are passages with which many readers will 
disagree. Let them state their: objections as clearly as Wiener stated his case, 
and they will have benefitted by their reading! 

Being a neurologist the reviewer is naturally inclined to evaluate the book 
peering out of his own little corner as it were. He derives some encouragement 
for his provincialism from his associations with Professor Wiener in the “Macy 
meetings” from which he carried with him the impression that the. author, too, 
is passionately interested in the brain and its working. The outline of a theo- 
retical neurology is there, but it is still so general as to require undoubtedly 
much further work before it will give us full insight into neurological problems. 
Such work has indeed been actively pursued by Wiener and Rosenblueth, and has 
even partly been published. Control of muscular activity by the brain is evi- 
dently achieved in different ways. The control of voluntary movement, which 
Wiener considers (see, especially p. 118) is achieved by a system of several feed-_ 
backs. The control of respiration, on the other hand, which Wyss and his col- 
laborators analyzed so brilliantly is obviously based on frequency filters in the 
solitary tract.: Hence, one may be well advised not to talk too glibly about “feed— 
backs.” , : 
How adequate Wiener’s definition of information will prove to be, only the 
future can tell. To define information as a simple choice between two alterna- 
tives should evoke some contradictions from the Gestalt psychologists. How far 
or when that symbolifie function which was emphasized by Craik (“thought 
models reality’) and which is increasingly commanding Attention in certain 
quarters, can be handled by Wiener’s theory, is another puzzling question. 

It is not the least value of a good book to set the reader thinking beyond the 
actual contents of-its printed pages. The preceding paragraphs intend first of 
all to show that Wiener did write a good book. In fact, any one who is or should 


GERHARDT VON BONIN 147 


be interested in problems of relations, be they of neurons or of human beings, 
or of any other units, or in problems of information, linguistic or otherwise, can 
do no better than to read faithfully not only Wiener but also all the references 
he gives! 

Unfortunately the publisher has done less than justice to the book. There 
are several misprints which any good proofreader should have caught. Thus on 
page 69 in formula 2.24 there appears “l.i.m.,”’ while formula 2.25 has correctly 
“lim.” On page 77, line 8 reads “note that is of the form and.” Obviously two 
formulae have dropped out. The lack of an index in a book which touches on so 
many subjects is also regrettable. But it should be stated in all fairness that the 
book is printed on quite heavy paper and seems to stand handling very well, 
hence it can be leafed through time and time again to find a lost fact or name. 

To have made Kybernetics (sic!) explicit and to have delineated its scope 
and its methods is Norbert Wiener’s great achievement. 


GERHARDT VON BONIN 
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